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HELICOPTER HIGH GAIN CONTROL 


Thomas B. Cunningham, Edwin C. Nunn 


SECTION 1 

INTRODUCTION AND OVERVIEW 


Introduction 


For most helicopters, some form of feedback stability augmentation is 
required to obtain acceptable flying qualities over the flight envelope. 

Recent designs include the use of model-following designs to meet certain 
mission requirements. Various command modes (such as attitude or velocity) 
can be designed with a model- following approach. Large feedback gains are 
desirable to improve tolerance of variations in mass and aerodynamic 
parameters, rejection of disturbances, and compliance with pilot commands. 
Whether these mo del- following designs are constructed using classical or 
modern control methodology, practical constraints limit feedback gain levels. 

Recognition of achievable gain levels must be part of the design process if 
high gain feedback is to be implemented. If gain limitations are known and 
are considered during the design, signal shaping and feedforward signals 
can be employed to compensate for some gain and bandwidth limitations. 
Neglecting these constraints during design will result in flight test problems 
forcing gain reductions which may seriously degrade overall system perform- 


ance. 



Benefits 


The primary benefit of high gain design is improved performance of VTOL 
aircraft through improved control. The prospect of reduced sensitivity to 
parameter variations is extremely important in the terminal area. 

Near-term benefits are: 

• Maximum bandwidth achievement of the three axis augmentation 
loops consistent with adequate stability margins 

• Identification of parameters limiting the bandwidth (such as 
actuator nonlinearity, rotor dynamics, etc. ) 

• Indication of ways of accommodating the gain limitations in 
control law synthesis to obtain satisfactory performance 

Overview 

The design study had two phases, modeling and design. Modeling was 
performed on the dynamics representing all gain- phase shaping elements 
of the vehicle and its control hardware. Hardware anomolies of gain 
limiting importance, such as control nonlinearities and sensor noise, were 
also examined. The control design effort was based upon successive 
eliminations of limiting factors imposed by modeled hardware constraints. 

In designs for the lateral- directional axes, comparisons to similar efforts 
performed at NASA Langley Research Center 1 were made. Finally, sensitiv- 
ity was examined at numerous portions of the flight envelope for necessary 
gain changing to achieve a consistent response in all regimes of flight. 
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Modeling 


The NASA /Army /Boeing Vertol CH-97B helicopter and its current control 
hardware mechanization were modeled in four areas. 

Rigid Body. - The aerodynamic data necessary to model the rigid body 

2 

motion of the CH-47B aircraft were supplied by NASA LRC. Results from a 
model verification flight test conducted as part of the study effort 
demonstrated close comparisons with the provided data. 

Actuator Dynamics . - Hangar and flight tests were conducted to identify 
actuator dynamics and nonlinearities and verify rigid body representations. 
Discussed in detail in Section 2 and Appendix C, these tests led to the 
development of dynamic models for the Electronic Control System (ECS) 
actuator and the lower and upper boost actuation systems. Also, important 
nonlinearities were examined and modeled. These included ECS rate limits 
and primary flight control system hysteresis. The pitch axis control 
hysteresis was extensive and eventually was the bandwidth limiting factor. 
Dynamic response impact was greatest from the upper boost roll-off 
characteristics. Second-order models for pitch, roll, and yaw upper boost 
actuators featured natural frequencies of 46, 50, and 37 rad/sec, respectively. 

Rotor Dynamics. - Flapping modes were modeled for each of the two 

rotors in tandem. Models were derived through multiblade coordinate 
3 4 

representation * and structured in state space to be appropriately mated to 

5 

the rigid body model. . Rotor interaction, an important consideration, was 

2 

modeled through empirically derived transcendental relationships. A 
twelfth-order dynamic model resulted, which allowed examination of rotor 
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feedback control. The model can be simplified to reduced order representa- 
tions. Coning for pitch axis control design reduced to a realization resembling 
a notch filter. Lateral and differential cyclic transfer through rotor dynamics 
reduced to fourth- order models. 

Sensor Models . - Rate and attitude gyro sensors were used for the rigid 
body feedback designs. Rate gyros displayed parasitic outputs at multiples 
of the rotor governor frequency (24 rad/sec). The magnitudes of these out- 
puts had to be reduced in all implementable feedback designs to avoid excessive 
actuator activity. Roll rate gyro noise, in particular, was high in magnitude. 
This proved to be the gain limiting factor in the roll control design. Verifi- 
cation and magnitudes for one- and three-per-rev noise models were derived 
from flight test results conducted as part of the contract effort and additional 
data supplied by NASA. 


Control Design 

High gain control allows design confidence in feed forward compensation to 
result in desired vehicle response characteristics. The specific goals used 
were high bandwidth second-order dominated responses for pitch and roll 
attitudes and first-order dominated response for yaw rate. 

Modal control, a new multiple input control design technique, 6, 7 was used to 
carry out the design goals. This technique is based on eigenvalue- 
eigenvector placement of the closed loop system. Therefore, the specific 
design goals can be embedded into the control design algorithm. 
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Pitch axis designs resulted in a closed loop control bandwidth which was 
limited by control hysteresis. Elimination of this will produce results only 
slightly better because of a pitch rate gyro attenuation filter needed to 
eliminate the one- and three-per-rev sensor noise. Eliminating the noise 
source or allowing higher actuator activity will allow higher bandwidth. 

This was examined by assuming no filtering on the pitch rate gyro. Rotor 
feedback allows higher bandwidth designs to be possible, but results in rigid 
body feedback gains too high to accommodate current sensor noise and 
actuator nonlinear characteristics. 

Roll and yaw control designs were conducted simultaneously with the modal 
control software. When using this approach, limitations were experienced 
only in the roll axis. Further expansion of yaw axis bandwidths for the cases 
studied was not tried. 

The fundamental limitation in the roll axis control is the one- and three-per- 
rev sensor noises on the roll rate gyro. Successful attenuation of these noises 
was achieved using notch filters at the dominant rotor frequencies, 24 and 72 
rad/sec. Although very limited band notches were used, the resulting phase 
losses resulted in little improvement over the complementary filter approach 
used at NASA.^ Further increases in the allowable bandwidth resulted from 
elimination of the gyro filtering. Rotor feedback again demonstrated benefits 
but again with extremely high rigid body feedback gains which would be 
impractical because of sensor noise. 

Gain scheduling to assure consistent response over the flight envelope will 
require simple airspeed schedules on series integral gains. 
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Throughout the design process, it was recognized that ECS rate limiting 
was a potential limitation. Evaluation of this would involve more detailed 
analysis of the feed forward input model. This effect was included in 
simulation models for transient response evaluation. Large input command 
steps resulted in system instability; therefore, rate limiting the feed 
forward commands would be required. 

Report Organization 


The report is organized to emphasize the most important results of the effort 
in the main body of the report, with detailed model analysis and testing and 
control design theory presented in the appendixes. 

Section 2 discusses the four modeling areas used: rigid body, actuator, 
rotor, and sensors. Section 3 contains the results for all control designs 
and evaluations. The effort is summarized and conclusions are drawn in 
Section 4. 

Appendix A contains modeling details for the tandem rotor flapping model. 
Appendix B discusses useful techniques for using modal control design 
theory. Appendix C contains details of the model verification flight test 
procedures and subsequent model construction. 



SECTION 2 


MODELING 

Introduction 

Feedback control design typically involves designing to uncertainties. The 
designer performs a tradeoff between closed loop performance versus 
sensitivity to model uncertainty. Performance can be either transient 
response, external disturbance rejection, or a combination of these. Model 
uncertainty can arise from a number of interacting sources: errors, 
variations, or omissions in the design models. A partial list for an aircraft 
would be: 

1. Rigid body model parameter uncertainty 

2. Known rigid body model parameter variations 

3 . Control hardware uncertainties 

• bandwidth 

• nonlinearities 

4. Unmodeled dynamics 

• flexure modes 

• flutter (more common to fixed wing aircraft) 

• rotor dynamics (for rotary aircraft) 
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5. Sensor constraints 


• bandwidth 

• parasitic anamolies, such as noise 

Rigid body model parameter uncertainty results from errors in accounting 
for aerodynamic characteristics. As a new vehicle design evolves, aero- 
dynamic data undergo continual upgrading from analytical models, through 
wind tunnel testing, to flight test results. A vehicle such as the CH-47 
helicopter has undergone periodic parameter identification throughout its 
lifetime. Uncertainty in its rigid body parameters, however, will never be 
zero, establishing one source of need for stability margins in the feedback 
control design. 

Known rigid body model parameter variations result from identifiable changes 
in a vehicle’s characteristics over its operating region. These characteristics 
can be complex nonlinear functions, but successful control designs have 
traditionally resulted using quasilinear approximations throughout a given 
vehicle’s flight envelope. Some adaptation is usually required to obtain 
desired performance over the range of known variations. Gain schedules on 
low frequency vehicle states are most common. Dynamic pressure, for 
example, is commonly used for fixed wing aircraft and forward speed for 
rotary wing vehicles. More elaborate schemes utilizing the basic flight 
control sensors have been successful?^*’ ^ 

Control hardware for aircraft usually refers to servo- actuator dynamics. 
Characteristics such as bandwidth and nonlinearities are usually known with 
greater precision than aerodynamic parameters. Hardware constraints are 
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the most common source of bandwidth limitation (academicians are notorious 
for ignoring these constraints). High gain control design in particular requires 
accurate modeling with respect to actuator dynamic characteristics. 

Unmodeled dynamics are referred to here as higher frequency vehicle 
characteristics. Flexure modes, flutter modes, and rotor dynamics are 
examples. Models exist for these dynamics but they can usually be charac- 
terized as 1) very complicated, 2) not very good, or 3) both. The easiest 
control design approach is to avoid these regions by low pass filtering, if 
possible, or by notch filtering to eliminate an uncertain band in frequency. 

Both techniques reduce the allowable bandwidth of the closed loop design. 

Active control approaches such as high frequency dynamic measurement and 
control are currently popular research topics. 

Rotor dynamic models for low advance ratio helicopters are in better form 
than either flexure or flutter representations. The reason for this is simply 
that rotors are more precisely fixed in operation by the vehicle designer. 

They are usually governed to operate at a constant rate throughout a vehicle’s 
flight envelope. This also affects greatly the vehicle's rigid body character- 
istics because low frequency aerodynamic characteristics are dominated by 
the rotor aerodynamics. 

Sensor modeling can usually be a secondary consideration because 1) the 
sensor bandwidth is typically very high relative to the desired design band- 
width (an order of magnitude in frequency is not unusual) and 2 ) wide and 
flat band low magnitude noise usually presents few problems. The sensors 
on the CH-47, however, contain very well defined high magnitude spikes at 
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multiples of the rotor frequency. These spikes produce one of the major 
obstacles to high gain control. 

CH-47B Rigid Body Model 


The NASA /Army /Boeing Vertol CH-47B helicopter is shown in Figure 1. 

Linear rigid body model parameters were supplied by NASA Langley Research 

Center. These were derived from nonlinear equations of motion which 

included vertical and rotor trim and rigid body state and control perturbations 

from the nominal trim values. Performed using a nonlinear model called 
2 

HELCOP, the perturbation-derived stability coefficients can be written in 
state space form: 


• A A 

x = Fx + Gu 
a a 


r = Hx + Du 
a a 


where 


x a is the aircraft rigid body state vector 


T 

x a = (u,v,w,p, q,r,0, 0) 


( 1 ) 


u is the control input vector 


u * <V V V V 


r is the response vector, in other words, typically all states and 

3 . 

controls plus other responses such as body accelerations 
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A 

F is the state compler matrix 

A 

G is the control input coupler 

A 

H is the response state coupler 

A 

D is the response control coupler 

Linear coefficient values for the rigid body are listed in Appendix E of 
Reference 2 for 99 trim points. The rigid body model is summarized in 
Figure 2. Some typical open loop performance numbers are shown in 
Figure 3 for various flight conditions. 

Rigid Body Model Validation 

As part of the study a flight test was conducted to determine actuator 
characteristics and validate rigid body model representations from HELCOF. 
Discussed in more detail in Appendix C, these tests included vehicle 
frequency responses at hover and 41.18 m/s (80 kn) forward speed. Selected 
results from the flight tests were compared with HE LCOP- derived frequency 
responses. These are shown in Figures 4-8. The agreement for data points 
less than 4 Hz is good. Beyond this point, however, the flight test data is 
corrupted by high magnitude sensor anomolies at multiples of rotor 
frequencies. 


Rotor Models 

The CH-47B rotor dynamics effect the high gain design limitations directly 
and indirectly. Direct impact results from the rotor physics producing 
attenuation. Indirectly, the impact is a result of noise spikes on sensors at 
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Figure 2. CH-47B Rigid Body Model, Taken from Reference 2, p. 89 
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multiples of the rotor frequency (24 rad/sec). This latter difficulty is 
discussed later. 


Rotor Flapping Equations of Motion 


Appendix A contains a complete description of the derivation of the coupled 
rigid body/rotor flapping dynamics for the CH-47B vehicle. Summarizing 
this development, the rigid body /rotor flapping dynamics can be modeled in 
state space 


x = Fx + Gu 
r = Hx + Du 


( 2 ) 


where: 

x is the rigid body /rotor state vector 



x is the aircraft state 
a 

T 

x = (u, v,w, p,q,r, 0, 0) 

a 


x is the rotor state 
r 

T 

X r " (B OF' B OF' B CF* P CF # B SF* B SF* P OR* P OR' P CR* 
B CR* P SR' P SR ) 


u is the control input vector [same as (1)] 
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r is the response vector 



r T - rotor states 
r 


F is the state coupler 



F is the state coupler from aircraft states to aircraft 

a* a , 

if 

equations 


F is the state coupler from rotor states to aircraft 
a, r 

equations 


F is the state coupler from aircraft states to rotor equations 

F is the state coupler from rotor states to rotor equations 

r» r 


G is the control coupler 



G is the control coupler to the aircraft equations 

a 

^ A A A A 

F Sj a t F and G a t G. F and G are derived by residualizing the rotor dynamics 

of equation {2 ) as explained in Appendix A . 
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G r is the control coupler to the rotor equations 

H and D are the response coupler matrixes from the state and control 
respectively. Equation (1) is augmented with rotor states to complete 
this. 


Transfer Function Rotor Model 

Although a twelfth- order model is necessary for individual rotor dynamic 
observance and (as discussed in Section 3) rotor feedback control, this model 
can be simplified without loss in validity for inclusion of rotor frequency 
impact on individual transfer functions. Figures 9 and 10 show the frequency 
responses of key transfer functions for the CH-47B. 

Pole- zero calculations yield the transfer function models shown in Figure 11. 
These can be appended to existing rigid body model transfer functions and will 
produce results very close to the state space models. Some observations can 
also be made about the flapping model used. 

• The collective model (blade coning) is a notch filter at the rotor 
frequency. This is supported by intuition, since coning has equal 
portions of up flap and down flap over one revolution of an input 
frequency at 24 rad/sec. The average value is therefore zero. 

• The models do not contain high frequency attenuation. 

• The transfer functions are almost invariant through the flight 
envelope. This is true despite significant changes in some rotor 
model parameters. The dynamics, however, are dominated by 
rotor frequency terms. 
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Figure 9. Analytical Model With Rotor Frequency Response (q/6 ) Hover 
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Actuator Models 


The major function of the model verification flight tests was to investigate 
actuator dynamic characteristics. Appendix C contains procedural details 
for this effort and a composite actuator model for the three input controls 
of interest. 

Figure 12 contains another version of the Appendix C model (Figure C-l). 
This version has been adjusted to mate with helicopter rigid body and rotor 
dynamics by establishing unity steady state gain, phase variable formulation, 
and a backlash model incorporated by means of an upper boost rate deadband 
term. This latter simplification is demonstrated in general in Figure 13. 

Sensor Models 


Sensor characteristics typically go unmodeled for flight control design 
because rough notions of bandwidth and noise are accounted for implicitly in 
the sensor choices. This is one of the few hardware areas where the control 
designer has some voice (because he is the only one with a need for flight 
control variable sensors). 

The flight control sensors for the CH-47B, however, display well defined 
vibration modes in multiples of rotor frequencies. In addition to the spectral 
plots displayed in Appendix C (Figures C11-C14) NASA has generated spectrum 
magnitude plots for sensor outputs. Shown in Figures 14, 15, and 16, these 
plots were generated without excitation of the vehicle. The 12 Hz spike in 
Figure 15 is particularly large. This occurs at three times the rotor frequency. 
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Figure 13. Hysteresis Model 

The actual source of this noise is probably related to rotor dynamics and, 
in particular, the fact that there are three blades per rotor. The rotor 
dynamics discussed earlier lead to an interesting conjecture about this. 

The notch filter description of the rotor transfer functions (Figure 11) leaves 
some doubt as to where the input energy at the rotor frequency has dissipated. 
The explanation of equal up and down flapping at the rotor frequency in the 
coning model also points to an off-axis response in cyclic, both longitudinal 
and lateral, for a coning command input. Since this phenomenon is not 
modeled, a matchup of this idea and the sensor output spikes could be explored. 

In any case, these oscillations corrupt the performance of feedback designs 
as gains on these sensors are increased. Specific details on how to deal with 
this problem are discussed in the next section. 
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SECTION 3 


CONTROL DESIGN 


Design Objectives 


The goals of this effort were to explore limitations to high gain control for 
the CH-47 vehicle and overcome these limitations where feasible within the 
constraints of the current hardware. The design is one which uses the 
current sensor complement and servo- actuation system. Additional designs 
which go beyond these constraints, such as rotor feedback and elimination of 
control nonlinearities and sensor noise, were explored for potential benefits. 

Benefits 

The primary benefit of high gain control is parameter insensitivity leading 

to simple control law modifications for numerous applications of model 

following. Figure 18 shows a simplified representation of this idea. If the 

controlled plant, represented by the multivariable series compensator H (s), 

s 

feedback compensator H (s), and vehicle, performs consistently throughout 

P 

the flight envelope, the feed forward model can be chosen with good predicta- 
bility of results. Also, if the closed loop response of the plant is "fast" 
relative to the model, then the feed forward model becomes a lag compensator 
in the sense that the model will be replacing fast dynamics with slower 
dynamics. 
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Figure 18. High Gain Model Following 

The objective of this study is to design the high gain feedback controls to 
enable simple model design for model following. Therefore, the high gain 
design does not concentrate on the model itself but attacks the problem of 
designing a fast consistent response throughout the flight envelope. 

Specific Design Goals 

The specific design goals are to construct high bandwidth second-order 
responses for 1) the pitch rate/pitch angle pair in the longitudinal axis to a 
differential collective input and 2) roll rate/roll angle pair in the roll axis 
to a gang lateral cyclic input. Also, a high bandwidth first-order yaw rate 
response in the yaw axis to a differential lateral cyclic is considered. These 
goals are shown in Figure 19. 
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Figure 19. Control Response Design Goals 
Modal Control Design Approach 


Modal control design techniques, outlined in detail in Appendix B, have been 
used here because of the consistency with which one can achieve design 
specificity in the face of 

• Higher order dynamics (rotors and actuators) 

• Coupled dynamics, for example, roll and yaw 

• Compensation, for example, proportional plus integral, 
notch filtering 

Also, the multiple designs necessary to explore the maximum bandwidths 
dictate the use of an inexpensive computation tool. The modal control design 
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software is very fast, thereby producing desired results with minimal 
computer expenses. Furthermore, the design procedure makes use of only 
specified outputs and requires no post-design gain adjustments. 


Compensation 


Helicopters typically have vibration components at multiples of rotor frequen- 
cies. For the CH-47 the most important characteristics are high magnitude 
disturbances on sensors, with narrow frequency bands centered at the rotor 
governor frequency and integer multiples of the rotor governor frequency. 
Such disturbances do not effect the controller response of the vehicle until 
high gains are introduced. Initially, the servo- system responds because of 
its high bandwidth. Eventually, the rigid vehicle response is effected by the 
rate saturation of the ECS servo, causing a net reduction in ECS bandwidth 
which in turn lowers the rigid body response damping. Two types of 
compensation can be used to handle this: 

• notch filtering 

• complementary filtering 

Notch Filters 

Because the undesirable noise occurs at such precise frequencies, a very 
narrow deep notch filter can be used to effectively attenuate the noise and 
preserve loop phase close to the notch frequency. 
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Complementary Filtering 


This concept has classically been utilized to combine sensors which contain 
information about a desired output in such a way as to attenuate individual 
sensor errors, including noise, and construct an output accurate over a 
broad range of frequency. 

The concept examined here is based upon earlier work done by Garren and 
Niessen. 1 A description of how a complementary filter is used in a feedback 
path is shown in Figure 20. 

Proportion Plus Integral Control 

Control design in all axes includes proportional plus integral compensation. 
The benefits of using this type of series compensation are 

• High gain is easily acquired at low frequencies. 

• The pitch axis loop design contains a low frequency zero 
which if left uncompensated dominates the transient response, 
producing high overshoot. This can be correct by either 
proportional plus integral control (as performed here) or a 
control input feed forward which cancels the zero in question. 

• Automatic trim adjustment results from proportional plus 
integral control. 
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WHERE 

He (S) = COMPENSATION FOR MEASUREMENTS y' m x m = MEASUREMENT OF OUTPUT, x, 

TO BE LOW PASSED 

y m = ALL MEASUREMENTS EXCEPT x m x = ESTIMATE OF OUTPUT, x, TO BE 

HIGH PASSED 

K x = GAIN FOR x FEEDBACK u = CONTROL INPUT VECTOR 


Figure 20. Complementary Filter Logic 




Feedback Control Design Verification 


Designs were verified using linear root locus techniques to establish 
dominant mode response for a given axis and assess the impact of all other 
system dynamics, including rotors, compensation, etc. Designs were also 
verified using step response simulation which included actuator nonlinearities, 
in other words, ECS rate limits and upper boost hysteresis. The simulation 
also included sensor noise spikes at one- and three-per-rev to test the 
effectiveness of notch and low-pass filtering. Results are discussed in the 
following sections. 


Fitch Axis Design Results 


High gain pitch axis control is achieved for various feedback arrangements 
for the differential collective. The following design concepts are based upon 
progressive constraint elimination. 

1. Pitch angle, 0, and low-passed (@ 2 5 rad/sec) pitch rate, q to 
differential collective, 6 . This case assumes current hardware 

-O 

complement, most important of which is the inclusion of rotor 
frequency pitch rate sensor noise. Also, upper boost hysteresis 
is examined. 

2. Pitch angle, 0, and nonfiltered pitch rate, q. Clean sensors and 
linear control hardware are assumed. 


3. Clean 0 and q plus front and rear rotor coning angle (3 and 

Or 

feedback. 


As discussed later, the ECS rate limit is a severe constraint; however, this 
is more application- oriented, since the demands imposed by the input model 
greatly affect the impact of ECS rate limitations. 
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Finally, the gain scheduling of design #1 is analyzed to produce a consistent 
control throughout the flight envelope. 

Pitch Angle — Low-Passed Pitch Rate Design 

Using the modal control software, the design for the pitch axis was developed 
by increasing the desired pole placement demand natural frequency, a> n , along 
a C = .707 radial in the S-plane. This is shown in Figure 21. 



Figure 21. Pitch Axis Design Goal 
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This is a single input problem, in other words, differential collective (6 ), 

33 

and therefore only pole placement was possible using the modal control 
algorithm. This produced no additional difficulties because the 9/6.^ 
transient response was dominated by the second- order root pair placed by 
the algorithm. 

The first-order low- pass filter at 25 rad/sec on the pitch rate sensor dis- 
tinguishes this design. The attenuation of the first and third rotor harmonic 
components of noise on the pitch rate sensor is thereby accomplished (see 
Figure 14). 

As with all design cases in this report, a root locus plot of design options 
best presents the design results in terms of essential performance. Results 
for the low-pass pitch rate feedback case are summarized in Figure 5. 
Included in the figure are the dominant roots of the system. The entire 
axis has 13 modes: four rigid body, four rotor modes (coning only), four 
actuator modes, and the first-order low-pass mode. In this case, however, 
only five modes are needed to examine effective performance because others 
are either higher in frequency or effectively cancelled by zeros. Designs for 
desired natural frequencies from 3 to 8 rad/sec are included in Figure 22. 

As the placed pole pair moves higher in frequency, a single pole from higher 
frequency dynamics moves into a position of dominance. The optimal set 
appears to be at a frequency of approximately 5 rad/sec. 


As shown in Appendix B, eigenvector placement requires more than one 
control input. 
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Transient response plots of two design cases tell a different story, resulting 
from the primary flight control system hysteresis (modeled as in Figure 12). 
Figures 23 and 24 contain differential collective step responses for two 
design cases. Figure 23 shows transient responses of various pitch axis 
dynamic states for the design case = 5 rad /sec and the effect of the 
hysteresis' on the performance. 

\ 

Appendix C outlines the derivation of the hysteresis model used. The pitch 
axis, shown in Figure C-3, is by far the worst axis, as is evident from the 
fact that hangar tests on this axis produced significant backlash. In general, 
there are a number of ways to attack a problem like this: 

1. Create a compensating nonlinearity in the loop to "cancel" 
the hysteresis. This approach is usually sensitive to changes 
in the fundamental nonlinearity, 

2. Provide a high frequency input dithering signal to excite the 
actuator in the dead region, thus naturalizing its effect. This 
approach is commonly used; however, in this case the ECS 
rate limit also constrains us to keep the activity low. Indeed, 
an ideal dither signal already exists in the high frequency rotor 
noise anamolies on the rate gyros. 

3. Redesign the control system to minimize hysteresis. 

One solution which helps is to reduce the gain on the pitch attitude. This is 
shown in Figure 24. Our transient bandwidth is reduced, but the impact of 
the limit cycle is much lower. This design is shown in Figure 22 as = 
7.0 rad/sec. In addition to the placed root pair at W D = 7.0, the response 
is also affected by the lower frequency roots shown. The net effect is a 
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40.000 50.000 






transient which displays characteristics of approximately 4.0 rad/sec in 
bandwidth and high damping. 

A stability check of this latter design shows 10.4 dB gain margin and 46° 
phase margin. We have, therefore, a good design for vehicle implementation. 

Pitch Axis Designs with Unfiltered Feedback 

As discussed in the previous section, the control system hysteresis produced 
the fundamental limitation to high gain design in the pitch axis. The next 
limitation was imposed by the sensor noise at rotor frequencies, in other 
words, the requirement of low-pass filtering. The next design sequence is 
based upon two assumptions: 

• Elimination of control nonlinearities, namely hysteresis 

• Clean sensors 

In Figure 2 5, the results of this design sequence are displayed. As the 
bandwidth demands increase, the system is increasingly dominated by a 
root excursion from higher frequency (originally part of a rotor pair). The 
design limit is = 8 rad/sec. This is compared to '=4.0 rad/sec 
for the hysteresis-dominated design and 5 rad/sec for the sensor noise- 
dominated design in which low-passed pitch rate was utilized. 

Rotor Feedback 

As demonstrated in Figure 2 5, the design limitation for high gains, assuming 
clean sensors and no hysteresis, was the rotor dynamics. Using the fourth- 
order coning model derived in Appendix A and the actuator dynamics, a rotor 
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Design Root Locus for Clean Sensors 







feedback design was analyzed. As depicted in Figure 26, only two modes 
are observable. This is also demonstrated with the second- order coning 
model presented in Figure 11. 


The feedback design of coning position states Pqjp and was based upon 
moving the lowest frequency roots of the rotor dynamics, located in W 


n 


24 rad/sec and C = .478, to a higher frequency. Within the linear constraints 

of the actuators, the lowest frequency attainable for the system of Figure 26 

was W =36 rad/sec with C = • 7. 
n 


The rigid body pole placement algorithm was then used as before except that 
the open loop plant contained the rotor feedback design. Figure 27 contains 
the results of this effort. As can be seen, the maximum achievable band- 
width is between Wj_ ) = 11 rad/sec and 12 rad/sec. 


Lateral- Directional Axes Design Results 


Roll and yaw control designs were analyzed using the coupled lateral- 
directional rigid body dynamics of the CH-47B. High gain control designs 
were emphasized for the roll axis; however, increased demands upon the 
yaw rate response were also imposed. As with the pitch axis, the lateral- 
directional axes design sequence proceeded along a similar scenario. 

1. Roll rate, p, yaw rate, r, and roll angle, 0, measurements 
were fed back to lateral cyclic, 6S, and differential cyclic, 

6R. Current actuation hardware, including rate limits and 
hysteresis models, is assumed. Rotor frequency sensor noise 
was included on roll rate and yaw rate outputs. 
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S-PLANE 



HOVER FLIGHT CONDITION 














2. Unfiltered sensors for p, r, 0. Clean sensors and linear 
hardware are assumed, 

3. p, r, 0, plus front and rear rotor lateral flap angles (p • 

SF 

and j3 ) are fed back. 

SR 

Finally, gain scheduling implications of design #1 are examined. 

Langley Complementary Filter 

As briefly outlined earlier, NASA Langley* has developed a unique approach 
to the high gain control problem. When sensor noise attenuation dominates 
the problem, as it does in the roll axis, the use of a "complementary" filter 
in the feedback loop yields good transient and noise rejection results. 

Figure 28 shows the specific arrangement of the NASA design outlined in 
general terms in Figure 20. 

This filter was examined here for its consistency of response for various 
choices of the high- and low-pass filter break frequency, W Figure 29 
shows a root locus of dominant poles for the roll axis for three choices of 
W c » These results verify earlier flight test results* which indicate similar 
loop responses for a similar range of W c values. 

A typical roll stick transient is shown in Figure 30. The low frequency 

A 

departure of the high frequency roll rate estimate, P, from the actual roll 
rate value, p, is indicative of the high-pass/low-pass nature of the design. 
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Filtered Roll Rate Design 


The problem attacked by the complementary filter is the same one presented 
in the current study, namely, attenuation of the one- and three- per-rev 
sensor noises. The roll rate gyro presents the biggest problem in this regard, 
because these noises, as shown in Figure 15, have the highest magnitudes of 
all the gyros (Figures 14 and 16 for pitch rate and yaw rate, respectively). 


The approach taken here is to attenuate these noises as much as possible 
and retain as broad a spectrum as possible of the sensor's output. The 
guideline for choosing the appropriate filter feedback combination is to limit 
the amount of sensor noise which is passed to the actuators. This limit was 
chosen to be less than or equal to the noise which passed through the 
complementary filter. The following filter provided this attenuation while 
maintaining a broad band of roll rate. 


TFC(s) = 


/ _2 a 2 1 
S +24 


1 S 2 + 72 2 ^ 


60 

\S 2 + 2(. 1) 24S + 24 2 j 


^ S 2 +2(. 1) 72S + 72 2 j 


S + 60 

l / 


The filter is composed of two narrow deep notch filters at one- and three- 
per-rev frequencies and a low-pass at 60 rad/sec. A 20 rad/sec low-pass 
was also applied to the yaw rate gyro. 


Modal control design was then applied to the filtered sensors. Because 
there are two controls, in this case, eigenvector placement can be performed. 
Figure 31 shows a selected set of eigenvectors for a typical design case. 

As outlined in Appendix B, the desired state responses are indicated by the 
"1.0" elements and the responses to be eliminated are indicated by the "0.0" 
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Figure 31. Eigenvalue- Eigenvector Placement Technique 



elements. The "-1.0" are coded by the algorithm to indicate arbitrary 
eigenvector mode responses. Pole placement is guaranteed, but eigenvector 
admissability is not. Comparing the left side of Figure 31 with the right 
indicates that our desires are met, on the whole. 

Figure 32 is a root locus of the eigensystem placement exercise for the 
filtered rate gyro cases. As with the pitch axis design sequences, the band- 
width expansion proceeds until a higher frequency root excursion moves into 
a position of dominance. In this case, this occurs between 4.0 and 5.0 rad/ 
sec natural frequency. A higher limit on yaw rate could be obtained but this 
was not explored further. 

The control law for the = 5 rad /sec case is shown in Figure 33. A 
lateral cyclic, 6s, step input transient is shown in Figure 34 to verify the 
second-order response of the roll attitude. Some hysteresis limit cycling 
occurs but this is much smaller than observed in the pitch axis, as expected. 
Also, a comparison of Figure 34 with the complementary filter transient of 
Figure 30 demonstrates only a slight improvement with the current design. 
Figure 3 5 demonstrates the first-order shape of the yaw rate response for 
a step differential cyclic, 6R, input. 

Lateral-Direction Design Without Filtering 

In order to gain some perspective as to the losses due to filtering and explore 
higher gains based upon more idealized hardware, that is, clean sensors, 
bandwidth expansion designs were conducted without sensor filters. 
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Figure 36. shows the results for this sequence. The extra sensor band 
allows the closed loop bandwidth to increase to a "p-0 " closed loop natural 
frequency of 7.0 rad/sec before troublesome root excursions from higher 
frequency start to dominate. 


Rotor Feedback 

Using rotor feedback to lateral cyclic demonstrated some additional band- 
width benefits. Using rotor dynamics derived in Appendix A with the lateral 
cyclic actuator model, a successful design was performed. As shown in 
Figure 37, the use of rotor feedback allowed the M P-0 " root pair to achieve 
a 10 rad /sec natural frequency. 

Finally, the gain values listed in Figure 37 are extremely high from a 
sensor noise standpoint. Since no post-design gain sensitivity was performed 
one can only conjecture that some may be arbitrarily high; however, most 
of those listed gains are much larger than those without rotor feedback, as 
was true for the pitch axis design with rotor feedback. 

Gain Scheduling 


The filtered designs for longitudinal and lateral directional axes were 
examined at forward speeds from -20. 5 to 61. 77 m/s (-40 to +120 kn) at 
zero vertical speed and -10. 16 m/s (-2000 ft/min) to + 10. 16 m/s (+2000 
ft/min) at zero forward speed. The modal control software produced very 
similar gains to the hover condition, thereby raising the possibility of a 
constant gain system. 
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Two difficulties with the transient responses surfaced: 

1. At extreme forward speeds (both negative and positive) the 
pitch axis response overshoot increased from approximately 
13% at hover to over 20%, This was corrected by increasing 
the integral gain in the pitch loop. A schedule on the integral 
gain with airspeed would be in order. 

2. At extreme negative vertical speeds the yaw rate response did 
not return to steady state within a reasonable time. Here the 
solution was to increase the yaw rate integral gain, thereby 
suggesting a schedule for this parameter. 



SECTION 4 


SUMMARY AND CONCLUSIONS 


\ 


The prospect of high gain benefits must always be tempered with hard 
reality. The approach taken was one of looking at a hierarchy of constraints 
beginning with those that exist, eliminating them one by one, and assessing 
the benefits. If a given limitation was solvable with the current hardware, 
the appropriate recommendation (such as compensation) was made. 

Other limitations involved characteristics or anomolies not likely to be 
removed easily. For example, sensor noise spikes are a fact of life for 
this set of hardware. The approach here was simply to legislate away the 
problem and examine the benefits of doing so. 

Modeling 


Modeling was a major program task. Rigid body data was supplied at the 

2 

outset in high quality form. Conclusions about other models developed on 
the contract are: 

• Actuator Models 

- The CH-47 boost actuators have considerable bandwidth, 
sufficient for all cases studied. Control system hysteresis 
posed difficulties in design performance, particularly the 
pitch design, where this nonlinearity ultimately became the 
bandwidth limiting factor. 
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- The ECS rate limit was not identified as a limiting factor 

in this study; however, this limit has the potential to impact 
the performance of all research system designs. This limit 
would greatly effect the response time if saturated. But more 
important, the system can go unstable if sufficient saturation 
is achieved. It is suggested that either this limit be increased 
or feed forward designs be constructed to avoid excessive ECS 
rate excitation. 

• Rotor Models 

- Good models for transfer function representation for first 
harmonic flapping modes were derived with minimum realiza- 
tion. A second-order "notch filter" type coning model and 
similar fourth -order representations for lateral cyclic were 
derived. These demonstrate little parameter variation through 
out the flight envelope. 

- Expansion of rotor models would logically include higher 
frequency flapping harmonic (three-per-rev) and lagging. 

This effect, however, should be preceded by more attention 

to unmodelled lower frequency dynamics such as the first body 
torsional mode. 

- Completely coupled state space representation of the first 
flapping mode was useful for rotor state feedback designs. 
Utilization of the reduced order transfer function realizations 
was not tried, but the use of similarly derived transfer 
relations from control inputs to flap states should also be 
successful. 
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Control Designs 


Control designs were based upon placing increased bandwidth demands upon 
a given set of assumed dynamics and finding the maximum possible band- 
width for a given axis (pitch and roll/yaw). Using second-order responses 
for pitch and roll as a guideline, the designs proceeded along damping 
ratio = .707 radials in the S-plane. A summary of the limits found is 
shown in Table 1. 

The roll/yaw axes design was compared with the NASA complementary 
filter 3 The conclusion here is that for the specific noise rejection problem 
which dominates the roll axis the complementary filter is a much simpler 
alternative and transient performance was equal to the current design. 

Finally, the use of modal control design techniques leads to some conclusions 
about this tool. Of the existing "modern" control design tools the one used 
here demonstrated the best facility for 

• producing practical, implementable designs without post- 
design gain changing and/or observer construction 

• conserving the computer budget by being computationally 
efficient 

• allowing one to input design goals that make sense; in other 
words, eigenvalue- eigenvector placement has direct relation- 
ships with classical design specifications. 
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TABLE 1. HIGH GAIN CONTROL"' '' DESIGN RESULTS 


Axis 

Controls 

Used 

Sensors 

Used 

Design 

Criteria 

Assumptions 

Filtering 

Bandwidth 

Limit 

Reason for 
Limit 

Pitch 

6 T3 

q, e 

2nd-order 

q-& 

response 

-Nonlinear 
actuators 
-Noisy sensors 

25 

flUer - S+25 

4, 0 rad/sec 

Upper boost 
hysteresis 

Pitch 

5 b 

q* e 

2nd-order 

q-e 

response 

- Linear 
actuators 
-Noisy sensors 

q filter g+25 

5. 0 rad/sec 

. Linear system 
performance 

Pitch 

6 b 

q> ^ 

2nd-order 

q-0 

response 

-Linear 
actuators 
-Clean sensors 

No filtering 

8, 0 rad/sec 

Linear system 
performance 

Pitch 

5 

B 

q, S 

9 of’ 9 or 

2nd-order 

q-6 

response 

- Linear 
actuators 
-Clean sensors 
-Rotor position 
states available 

No filtering 

11.0 rad/sec 

Linear system 
performance 

Poll 
k Yaw 

6 6 
S’ R 

p» r ( 0 

-2ncI-order 
p-0 response 
- Ist-order 
r Response 

-Nonlinear 
actuators 
-Noisy sensors 

p filter * 

S+24 2 S+72 2 

5.0 rad/sec 

Sensor noise 




S 2 +4, 8S+24 S 2 +14, 4S+72 2 








60 

X S+60 








20 

r filter = 

S+20 



Roll 

k Yaw 

5 , 6 
S R 

p t r, 0 

-2nd-order 
p-0 response 
-Ist-order 
r Response 

-Nonlinear 
actuators 
-Clean sensors 

No filtering 

7.0 rad/sec 

Linear system 
performance 

Roll 

k Yaw 

5 . 

S R 

p, r, 0 

a p 

SF' SR 

~2nd~order 
p-0 response 
-Ist-order 
r Response 

-Linear 
actuators 
-Clean sensors 

No filtering 

10.0 rad/sec 

Linear system 
performance 


Limit based upon roll performance. Yaw bandwidths for each design was a first order root at the indicated location 
Proportional plus integral control was used on all designs 




APPENDIX A 


ROTOR DYNAMICS 


Introduction 


A significant portion of the study effort involved modeling, that is, vehicle 

rigid body, actuator dynamics and nonlinearities, sensor noise, and rotor 

dynamics. Rotor flapping dynamics were modeled over one period of 

3 

revolution, as outlined by Hohenemser. Equations were patterned after a 

4 

development by Hall. Significant impact from the tandem rotor arrange- 
ment (in particular, the front/rear rotor interaction) required special 
treatment as discussed in reference model parameters provided by Ostroff, 
et al? 

The final model derivations are part analytical (kinematics and coordinate 

rotations) and part numerical (external forces, rotor interaction, and rigid 

body /rotor coupling). The rotor models are mated to linearized rigid body 

dynamics with a reverse residualization technique to make full use of 

2 

previously derived rigid body data. 

Blade Equations of Motion 


A three degree of freedom (flap/lag /pitch) dynamic model of a helicopter 
blade is developed in this section. 
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Coordinates and Rotations 


The coordinate systems used in the derivation are pictured in Figures A-l 
and A-2 and defined in Table A-l. The fuselage coordinate system (F) is 
fixed in the helicopter, centered at the vehicle center of gravity. The hub 
system (H) for each rotor is fixed in the vehicle, but centered at the respec- 
tive rotor hub and aligned with the rotor shaft which is tilted about the Y 

X? 

axis through an angle -i as shown in Figure A-l. The shaft coordinate sys- 
tem (S) is a fixed rotation (+180° about y ) to aid in completing the rotations 

IX 

necessary to derive equations in the blade system (B). The other axes sets 
all rotate with the blades. All orientations with respect to each other are 
shown in Figure A-2 and defined in Table A-l. 

Force and Moment Equations 

The equations of motion for the blades are derived from the following momen- 
tum relation 

£M = ft + x H + mR xa° _I <A-1) 

cm 

where, 2 M is the vector of external moments acting at the blade hinge point 

H is the angular momentum of the blade 
-> 13-1 

Q is the angular momentum of the blade 

m is the mass of the blade 
* 

R is a vector pointing from the hinge point to the blade center of mass 
cm 
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R = 9.144m 











TABLE A-l. ROTATION SEQUENCE 
(Fuselage to Blade) 


Rotation 

Axis of 
Rotation 

\ Axis System 

Frame Symbol 

Angle J " 

F 

(Fuselage) 

Tilt 

y f 

\ , 1 , ZL 

• \. JL JL JL 

/X Y 7 

H 

(Hub) 

Rotor Align 

y h 

“aioA H ' H ' H 

/(X„. Y„, ZJ 

S 

(Shaft) 

Azimuth 

Z 0 

s 

, \ kJ K-J K-/ 

* \ 

/(X Y . Z ) 

R 

(Rotor) 

Flapping 

y r 

^ R R R 

/(X tf Y_, ZJ 


Inplane Lag 

Z B 

- c r b b 

/(X r , Y r , ZJ 

c 

(Inplane lag) 

Pitch 

X 

e \ C C C 

/(X D , Y ZJ 

B 

(Blade) 

— — ^ U p u 


a°~* is the acceleration of the hinge point with respect to the inertial 
coordinate system 


The angular momentum of the blade is the product of the inertia matrix and 
the angular velocity vector. 


H 


-»B-I 
= i n 


(A-2) 


where 



I 

XX 

I 

xy 

0 

t— 1 

HO 

I 

xy 

I 

yy 

0 


0 

0 

I 




zz 


74 



I and I inertia products in this matrix are assumed zero due to the 
xz yz 

smallness of the blade thickness relative to the chord and span. 


”*B”I 

The angular velocity vector Q is given by 


n B_1 = n H_I + (?*' H + n B ' R 


(A-3) 


where. 


->H~I 

Q 


is the angular velocity of the hub with respect to the inertial 
axes. Written in hub coordinates it is 


n H_I - r 
Q = [p. 

H 


H’ q H’ r H ] 


Q is the angular velocity of the rotor with respect to the hub. 
Expressed in rotor coordinates it is 

4T) „ TJ rp 

n = [o, o, n] 

R 


B *R 

Q is the angular velocity of the blade with respect to the rotor 
coordinate system. It is a combination of blade flapping, 
lagging, and pitching motion. In blade coordinates it is 


^B-R 

Q 

B 



[9, 0 , 0] T +T 0/C 


T c/e (0 ' ' C > T 


+ T e/C ’ T C/r T B/R [0 * 4 ‘ 0]T 


The T y. matrices in this expression are transformation matrices from 
the j coordinate system to the ith coordinate system. 
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If all the terms in Equation A-3 are expressed in the blade coordinate system 
and assuming the flap, lag, and pitch rotation angles (0, C» 0) are all small, 
the result is 


->B-I 

o 

B 


-p cost + q sin'll + 00 + 0 - CP 
hi hi 

p sin'll + q„cost + 00 + 0£ - 0 - 00 
H H 

-r + Q + C + 0 0 
hi 


(A-4) 


Thus the angular momentum of the blade expressed in blade coordinates is 


I (-p cos'll + q„sint + 00 + 0 - CP) 

XX hi hi 


H-in 8 ' 1 

B B 


+ ^(p^.sin'li + q^cost +0O + 0C - 0 - C P^) 

I (-Pacos'!! + q sint + 00 + 0 - CP) I (A-5) 

xy rl Jtl 

+ I (p TT sin'li +q TT cos'li + 00+ 0*C - 0 - CP n ) 

yy h h 

W- r H + n + t + 0 ® ) 


The acceleration of the hinge point a° * in Equation A-l is given by the 
equation 


-»o-I ->h-l *R-I ->o-h *?R-I ,nR-I ->o-h. 

a =a +0 xp +0 x(0 xp ) 


(A-6) 


->h-I 

where a is the acceleration of the rotor hub 


0 


R-I 


is the angular velocity of the hinge with respect to the inertial 
axes. In rotor coordinates. 
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Q 

R 


R-I 


-P 


H 




cos^ + q sin\|< 
H 

sinty + q cosilr 

w 

+ n 


-»o-H 

P 


is a vector pointing from the rotor hub to the blade hinge point. 
Expressed in rotor coordinates it is 


H = [e, 0, 0] T (See Figure A-2) 

R 

Similarly a , the acceleration of the hub, is given by 


U H +q HV r H V H 

V H + r H U H “ P H W H 


-»H-I 

a 

H 


W H + P H V H " q H U H 


The acceleration of the hinge point a° * is found by solving Equation A- 6, 

taking care to transform all terms into the same coordinate system. If all 

X 

the small terms are dropped, the resulting expression for a is 
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-(u H + q H w H - r H V H ) cos ' t ’ 

+ (v H + ■ P H W H ) Sin '' 1 

o 

- e(n + 20r ) 

H 


-> o-I 

a 

B 


(1 H + q H W H • Vh 1 s 10 * 

+ + r H U H " P H W H ) COS * 

2 2 

- er + e£ 0 + eg 2Q 

H 


" (W H + P H V H " q H U H ) 

- e(p H - 2q H Q) sim|r - e (q H + 2p H Q) cos<|i 

+ peQ 2 - C9en 2 


(A-8) 


The vector R is Equation A-l. Written in blade coordinates, this is 
cm 




R 

B 


cm 


[6, T|. 0] T 


(See Figure A-2) 


(A-9) 


Substituting Equations A-4, A-5, A-8, and A-9 into Equation A-l and 
performing the indicated operation results in the general flap-lag-pitch 
equations of motion. 
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Pitch: 


M = I [{-p *q 0) cos ^ + (q + p ft) sint + pQ + 0 -CP] 
x XX H H K H (A-10a) 

+ (I - I ) [p„n sinijr + q Q cost + 9fi 2 + 20Cfi - £n - CP^ 2 ] 
zz yy 1 H H 

+ I [(P H " 2q H Q) sin t + (q H + 2p H Q) cost + e’£ - P - PQ 2 - 2p£0] 

+ m [-(w + p v - q u ) - e(p - 2q Q) sint 

i| n h n n H ii rl 

2 2 

-e(q + 2p 0) cost + PeQ - C e eQ ] 

xi ri 

Flap: 

= lyy[(P H - q H n) sint + (q H + P H n) cost + 9Q + e£ - P - CPO - Pen] 

+ (I - I ) [-p Qcost + q Q Sint + PO 2 + %Q + BQ - £*pO] (A-lOb) 

xx z z tl il 

+ I X y[("P H +2q^Q) cost +(q R +2p H Q) sint +9 - CP+ 9^ 2 +20£Q- CP^ 2 ] 

- mi[-(w H + P H V H " q H U H } " e( PfT 2 q H n) sin * ' e(q H + 2 p n Q) cos * 

+ peQ 2 - C6eQ 2 ] 
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Lag: 


M = I [-r + C + 90] 

z zz 1 H 


+ (I - i ) [-PP + een] 

yy xx L 

+ I [2 pen + 2 e'en] 

xy 


(A-lOc) 


+ m6[(^ + q H w H ’ r H v H ) cos ^ + ( V H + r n U H " COS ' 1, 


• 2 2 , 

er + + e0 0Q ] 

xi 


- rn [-(u H + q R w H - r H v H > cos t + < V H + r H u H - P H W H ) sini[r 


- eQ + 2eOr ] 
ri 


Some simplification in these equations is possible with the assumptions that 

I << I ~ I 
xx yy - zz 

I - I ~ I 
zz yy - xx 

The first of these assumptions implies that the moment of inertia about the 
flap hinge is approximately equal to the moment of inertia about the lag 
hinge and much larger than the moment of inertia about the blade pitch axis. 
The second approximation comes from assuming that the blade thickness 
and chord length are small in comparison with the span. 

Furthermore, if it is assumed that the blade is a homogeneous rectangular 
prism, its inertia properties can be approximated as 
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i ~i - m(R : e)2 
zz - yy 3 

r _ m(R-e)Tl 

xy " 2 

where m is the mass of the blade 
R is the blade radius 
e is the hinge offset distance 

T] is the distance from the pitch axis to the blade center of gravity 
(see Figure A-2). 

Substituting these relations into Equations A-10 and introducing the notation 

Ti = Tl/R 
e = e/R 

gives the simplified flap-lag-pitch equations of motion. 


Pitch: 

M 

= (~P H + 2q H Q) cos* + (q H + 2p n n) sin* + 0 + (Q +2Q £ ) -£(& +Cl$) 
xx 

I 

+y^- [(P H - 2q H Q) sin* + (q H + 2p H Q) cos* + £8 - 

XX 

m 

+— l-(w H + P h v h - q H u H ) - e(p H - 2 q[j a) sin* - 
XX. 

+ eQ 2 8 - eQ 2 £0 ] 


P -(n 2 + 2QC)e] 

(A-lla) 
e(q R + 2p H Q) cos * 
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Flap: 


r-i — = (Pu- - 2q Q) sint + <q + 2p Q) cos t - 1 - (Q + 2fi £)0 + 6 £ 

I Jtl xl no. 

yy 

+ J (1 _^ [{-P H + 2q R 0) cost + (q H + 2p H n) sint + 0 (A-llb) 


+ 


+• 


+ (0 2fi£)e - (0 + 0 P)£] 


3/R , 

2(1 --5) fW H P H V H “ q H U H ] 

» 2 2 

2(1 .^ [(P H - 2q R n) sint + <q H + 2p H Q) cost - n B + Q £ 0 ] 


Lag: 

M 2 

—■£. = -r + £ + el - pen + een + cr p e + eg 

I H 

zz 

+ wW tznee +20 66 ] 

3 / R • * 

+ [(U H + q H "H ' r H V H ) 9int + (V H + r H U H ' P hV C ° S * 1 

+ — li-g- [(U + q 0) - rv) cost - <v + r u - p w ) sint ] 

lx XX *0» XX JljL JTX XX XX XX XX 

-e) 

+ W^T t-'H + sn2+ pen2] + jrh - 2 nr H> (A - Uc) 

\ x ** 6 / 
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External Moments 


The external moments acting on the rotor blades come from two sources, 

(a) aerodynamic forces and (b) spring forces in the hinges. The moment due 
to aerodynamic forces is found by integrating the lift force acting on the 
blade. Lift per unit span is given by the equation 

1 2 
L, = paeU a 

where L is the lift force per unit span 
p is air density 

a is the slope of the sectional lift curve 
c is chord length 

U is the velocity of air passing over the blade 
a is the angle of attack of the blade 

If the radial component is ignored, the velocity of the airflow at the blade 
elements can be resolved into a component in the plane of rotation, U T , 
and a component perpendicular to the plane of rotation, U_, as shown in 

sr 

Figure A-3. Since U /V_« 1, the total lift acting in the vertical direction is 

P ^ 

L rr ac / (u T 2 + u p 2 ) ( e-u p /u T )dr 

R 

-“-g-pac J* U t 2 (1 + U p 2 /U T 2 ) (6 - U p /U T ) dr 

~ 1 n ^ 2 

- ~2 pac J <U T e - U T U p ) dr. 

0 
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4 

Figure A-3. Resolution of Blade Element Velocity Components 

The moment acting about the flap hinge line (-Y direction) is 

s 

~ 1 ^ 2 

M p = y pac J* (U T 9 - U T U p ) r dr <A-12a> 

and the moment about the lag hinge line (Z^ direction) is 
1 R 2 

M c = y pac J* (U p - U T U p 9) rdr (A-12b) 

where 0 is the pitch angle of the blade. 

With the introduction of the rotor lock number 
4 

= pac R 
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where I. is the moment of inertia about the flap hinge (T = I = I ) becomes 
p p yy zz 


M 

y 

I 

yy 

M.. 

z 

I 

zz 

where 


= Y^ = T J* 4(U^ 9 - u^u^xdx 


0 


P T 


- M c _ Y [‘ 

-r O J 


4(U p - UpUp 6) xdx 


(A-13a) 


(A- 13b) 


u T = u t /r 
Up * Up/R 


r/R 


Expressions for Up and Up are obtained from the equation 

-> B-I ->S - 1 -? R - 1 -> -*B-I -»b-0 

v = v +0 xe+Q x 6 

-»B-I 

where v is the velocity of the blade. In the blade coordinate system it is 

v B_I = [U R , U T , U p ] T . 

S- 1 

v is the velocity of the shaft. Expressed in hub coordinates it is 

-»S-I . ,T 

V “ [U H' V H W H * 

H 

R _ I 

is the angular velocity of the rotor with respect to the inertial 
coordinate system 
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p__ cost + q sin^ 

hi ti 

P H sint + q R cost 

r H +t2 

e is the vector from the rotor hub to the hinge line 

e = [e, 0, 0] T . 

R 

I 

^ is the angular velocity of the blade (see Equation A- 4). 
y b ~ 0 

o is the vector pointing from the hinge point to a point located a 
distance r down the blade 

t b “° = [r, 0, 0] T . 

B 

Performing the algebraic manipulations indicated in Equation A- 14 results 
in the following expressions for Up and Up: 

U T = (r+e) (n-r R ) + rC + (v R + Cu^.) cost 

- (u R + Cv R ) sint (A-15) 

U p = rB + [-(r+e) q R + pu R ] cost + [-(r+-e) p R - &v R ] sint - w R 

The total velocity at the blade is obtained by adding a downwash term to 
Equation A-15. Analytical and empirical formulae exist for downwash and 
these may be sufficient for single main rotor vehicles. However, for a 
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tandem rotor helicopter with overlapping horizontal rotor planes the down- 
wash velocity is the source of key cross rotor interaction terms. Rotor 
interaction will be discussed in a later section and its inclusion here is 
limited to the following definition 

v is the blade downwash (or induced) velocity generated by the 
increased air passing through the rotor producing thrust 

v is therefore added to the Up equation in A-15. 

U T = (r+e) + r£ + (v H + £u R ) cost 

- (u H + Cv R ) sirnlr (A- 16) 

Up = rp + [-(r+e) q H + PUp] cost + [-(r+e) pp - pVp] sint 

- w + v 
ri 

It is further assumed that the pitching motions of the blade are small, so 
that for a given blade element, the pitch angle of the blade, 0, can be 
approximated as 

9 = 9-0 sint - 9 cost + 9-r +h.o. (nt)' (A-17) 

US C 1 

where 9^ is the collective input 


Higher order harmonics of t can be used if appropriate. Hohenemser 
has expanded these terms inRef. 3. 
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°s is the longitudinal cyclic 'input 

•sV 

8 c is the lateral cyclic input 

0 T is the blade twist per radial element (assumed constant) 

This approximation also implies that the actual pitch angle dynamics 
(Equation 11a) are fast and therefore will be ignored. 


Spring Moments 

Spring restraints on the rotor flapping and lagging can be expressed as 
Mg spring = K^p 

(A-18) 

spring = K^6 

Linear Solution to Tip Path Plane Equations 


As with the pitch angle, the use of the first harmonic approximation is 
applied to the flap equation A- lib. From this point on, pitching, ©, will be 
assumed to be an input (Equation A-17) and lag, 6, will be assumed negligible, 
that is, 6 = 6 = 5 = 0. 


Letting 


P(t) = P (t) - P (t) cos^ - P sinf 
o c s 


(A-19) 


Special note should be taken of the terminology of longitudinal cyclic reference 
to the sin-| coefficient, 0 C . Control inputs effect blade flapping gyroscopically; 
that is, a sin'll direction input produces a cos'll output in 0; likewise, a cos'll 
input 0 produces a sin^ output in P. 
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where 


0 (t) is the coning angle 
o 

P c (t) is the longitudinal flap angle 

0 (t) is the lateral flap angle 
s 

and taking derivatives of A- 19 yields 

0<t) = 0 - (B„ + 00 ) cos* - (0 - 00 ) sin* (A-20) 

o c s sc 

B(t) = *0 - <0„ + 2 00 - Q 2 0 ) cos* (A-21) 

U C S C 

.. • 2 
- (0 -2 00 - 00) sin*, 

s c s 

Substituting A-19, A-20, and A-21 into the flapping Equations A-llb, we 
get grouped terms of rotor harmonic coefficients as follows: 


Coning: 

f 3“ g K q 3/R 

+ | [1 n + y J B 0 "2 (TIT <"h + Wh + u ho> 

n _ (A-22 ) 

= T M .6x + M (x ) 
oi l oo o 

i=l 

Longitudinal flapping, that is, sin(*) terms: 


** 2 i 3e 

0 +200 - 0 0 +< [1 + 
c s c ' 1 


2 K B 

2(1-1) ln *T- K 

yy 


Q p, Q p i 

+2n + 2 <r^r ]q H 


(A-23) 


n 


= T M .6x. + M (x ) 
ci i co o 


i=l 
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Lateral flapping, that is, sin(i|f) terms: 


9 s- 2nB c- n2 f S + ([ 1+ 2^)l n2 + 


yy J 


Q p Q p • 

20 [1 + 2(Cl-e) ] % + [1 + 2(l-e) ] P H 



M ,6x. + M (x ) 
SI 1 so o 


(A-24) 


where 

M (x.), M (x.), and M (x.) are trim moments for flap 
dynamic coning, pitching, and rolling, respectively. 

x q is the trim state vector containing trim values for all rigid 
body and rotor states and controls. These values are derived later. 


M q ., M c -, and are aerodynamic moment perturbations for 
state x. for flap coning, pitching, and rolling, respectively. They 
are derived below. 


5x. is the perturbation of state X. when X = (u , v , w , p , 
1 # # 1 H rl H H 

q^., r , P , P , P , p , p , P , 0 , 9 , 0 ) 
HHooccssocs 


Aerodynamic Moment Perturbations 


Analytical expressions for the individual coefficients for the M .'s, M .‘s, 

3 4 l ci 

and M . have been derived by Hohenemser and Hall. For this study, 
si 

however, these parameters were derived numerically because of numerous 
trim values. This created a massive set of perturbed combinations and a 
complicated, yet important interaction between the front and rear rotors 
(detailed later). 


90 



Using Equation A- 13a, 


M ^ 1 

y g v 

I 8 

yy o 


/ 


4(U t e - u p u T ) xdx 


and neglecting the offset distance, e. 

Up = (0 - r H > + (v^cost + Up.siniji)/R 

u pv 

Up = xB - (xq H - -J^-) cos if - (xp H +— ^-) sin^i 


- (Wp + v/R) 


(A-2 5) 


A- 13a can be expanded analytically, and certain major portions of the 
integral have been derived through first and higher harmonics in ^ . For the 
current application, however, the integration with respect to x is carried 
out by grouping terms in powers of x and performing the integration analyti- 
cally. This leaves numerous harmonic terms associated with 8 (Equation 

* 

A-17), 8 (Equation A-19), and P (Equation A-20), plus the harmonic terms 
contained in the definition of Up and Up (Equation A-25). The coefficient 
grouping is performed by numerical harmonic matching; in other words, 

2tt 


M 


M 


M 


2 TT 


2rr 


1 

2tt 


I 


0 


I 

0 

I 


2tt 


2rr 


M 

y 

l yy 

d^ 

(A-26a) 

M 



y 

r yy 

cos^d^ 

(A-26b) 

M 

y 

i 

yy 

sin(id^ 

(A-26c) 
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Trim values M , M , and M are calculated by fixing the input state x 
oo co so J 6 ^ 

at trim (x ). Perturbation quantities, M . 's, M .'s, and M .'s, are 
o n 01 ci si 

calculated by forward and backward perturbing of the input state. For 
example 


1 f 2tt M 1 y* ^ TT M 

2 -r, J I 2 - (x i + 4x i> cos * d * -2 ; J r 2 - (x i - ix i' 

J yy J 


2tt M 

yy 


cos'jid't 


M . = 

ci 


2 Ax. 


These results are coupled with the analytical relationships of Equations 
A-22 through A-24 to form the flapping equations of motion. 

Rotor Forces and Moments on the Vehicle 


The rotor impact on the vehicle comes from three sources: 

• Spring forces 

• Inertia forces 

• Aerodynamic forces 


Spring Forces 

Equation A- 18 shows the spring moment equations. Using only A- 18a, 
since we have dropped the rotor lagging dynamics, we can expand the flapping 
moment using Equation A- 19 for 8: 

MPspring = K n (P - P cos^ - p sin'll) 

P o c s 
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In hub coordinates this moment becomes 


M 


XSH 


M 


YSH 


= -MBspring sinty 
= -MBspring cos f 


Averaging over one revolution and including all blades (N ), 

jB 


N 


M. 


B 


XSH 2tt 


I 


2tt 


N, 


m sxh 


K-B c 
B S 


N. 


Mysn 


B 


/ 

0 


2tt N r 

M d* = _JL k B 
SYH v 2 Be 


Inertia Forces 


(A-27a) 

(A-27b) 


Using techniques documented earlier, the one- revolution average forces 
produced by inertia at the hub can be written in hub coordinates 

F X ih * -N b P b R [u h +q H (0 H - r H U H ] (A-28a) 

^T3 F [^tt Yj'^'LT ” Ptt^tt ] (A-28b) 

i Irl id jd il n. ix n il 

F ZIH “ ~ N B P B R ^ W H + P H V H " q H U H ~~2 " g o R ^ (A-28c) 

Moments acting on the fuselage at the hub result in the above forces being 
applied through the rotor offset e. These are derived in rotor coordinates, 
converted to the hub system, and averaged over one revolution to become 

Mx IH = " T eN B P B R [( ’^s +P H } ' 2Q( h +q H ) (A-29a) 
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(A-29b) 


M.„. rT . = 

1 

e N 

pR 

YIH 

4 

B 

B 

M„ = 

1 

e N 

p b r 

ZIH 

2 

B 


*♦ * 


(A-29c) 


The above equations also include the assumptions that products of small 
perturbations can be eliminated and e«R/2. 

Aerodynamic Lift Forces 

Equation A- 12 shows the lift characteristics for a given rotor blade. This 
can be translated to hub coordinate forces and moments as follows: 


F ZLH = " l 

m xlh “ e l sin * 

M YLH = " eLC ° S * 

where Equation A- 12 is modified to be: 

1 


(A-30a) 

(A-30b) 

(A-30c) 


_ 1 3 

L = y pac R 


/ 


(U T 6 - U p U T ) dx 


(A-31) 


0 


The A-30 equations are averaged over one revolution for all blades of a rotor; 



2rr 

f 2rr 



F = - 

ZLH 

J 

0 „ 

L d f 

(A-32a) 


“xlh" 

N e 
B 

- 2tt 

L sintydiji 

(A-32b) 

2tt 

J 

0 „ 





N p e 

2tt 

r 2rr 



^LH 


L cos^df 

(A-32c) 


0 
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Inplane Forces 


Figure A-4 augments Figure A-3 by showing the inplane force relationships 
for a typical rotor. In hub coordinates these can be resolved to be 


F = -D sin* + L(B cost - U_JU_ sin*) (A-33a) 

AJJ sr 1 

F yD = -D cos t - L (U p /U T cost + sin*) (A-33b) 


where 

D 


L is defined in Equation A-31 and 
1 3 r 1 ~ 2 

= T pCR C do J U T dX ‘ 


0 


(A-34) 


These forces are integrated over one revolution to yield the average value 


for an N„ bladded rotor, 
B 



F A f 

XDH 2n J 

- 2tt 


(A- 3 5a) 

0 

F [ 

YDH 2tt J 

■ 2tt 

d* 

YDH 

(A-3 5b) 


0 


Total Forces and Moments 


The forces and moments derived from spring, inertial, and aerodynamic 
sources are 


F 

= F 

+ F 

XH 

XIH 

XDH 

F 

= F 

+ F 

YH 

YIH 

YDH 

F 

= F 

+ F 

ZH 

ZIH 

ZLH 


(A-36) 
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“XH = “XSH + “XLH 


M = M + M 
YH YSH YLH 


m zh = 0 


(A-37) 


Transfer from Hub to Fuselage Coordinates 


In all previously defined equations the center of action was the rotor hub. 

All input dynamic states from the vehicle (vehicle translational and rotational 
velocities and accelerations) were input at the hub position. All output rotor 
forces and moments were also centered at the hub. 


It is necessary to translate and rotate these quantities to the vehicle body 
axis located at the c.g. of the aircraft. Referring to Figure A-l for a given 
rotor inclination, i, forward offset, l, and vehicle offset, h, the following 
vehicle state transformations of velocities must be made: 



(A-38) 


(A-39) 
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Since i, i, and h are fixed in time acceleration states, u H , v^., w^, p^., 

q , and r would be transformed in the same manner. The output forces 
H H 

and moments are transformed back to the vertical body system at the c.g. 
as follows: 



Rotor Interference 


The treatment of rotor dynamics, to this point, has proceeded nicely with 
individual rotors. The impact of rotor interference, however, is significant 
for this type of tandem arrangement. A good analytical study of this effect 

for the CH-46 helicopter is contained in Reference 5. 

\ 

The parameter v in Equation A- 15 represents the appropriate location for 
analyzing this effect. 

Let 

X T “ I — j — (interference in flow ratio) 

I ]U X\ 
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also 


w. 


X = 


M - 


H 

[qTr 


u. 


H 
[Q | R 


(in flow ratio) 
(advance ratio) 


Reference 2 contains the following front (subscript F) and rear (subscript R) 
interaction relationships 


'TF 


IF 


2 [(Xp “ ^ IF^ ^F ^ ^ 


2~ ! .2 ,1/2 + dF RF 


"TR 


(A-42) 


2 [<^ ' x ra )2 + ^r > 2]1/2 


IR 


2 t (X R " X IR ) (m R ) ^ 


TR +dF 

~ 2 ~. *2,1/2 dF FR 


'TF 


(A- 43) 


2[(X F - X lF ) 2 + (m p ) 2 J 1/2 


where interaction parameters dF „ and dF FD are defined empirically. 

Kx 1 Jt* R 


> 

For U = 0: 


dF pR = [0.356 + 0.321 R if - 0.368(R IF ) 2 + 0. 392(R IF ) 3 ] (1 - | sin P^,|) 


IF 


IF' 


(A-44) 


+ [0.356 + 0.0131 R t _ - 0. 0764(R ) - 0.0085(R X _) ] | sin pi, I 

I Jt* lx* I b r 


dF RF = [0.356 - 0. 151 R ir - 0.314(R IR ) 2 + 0. 164(R IR ) 3 ] |sin p^ | ) 


IR 


(A-45) 


+ [0.356 + 0. 0131 R - 0.0764(R__) 2 - 0. 0085(R T _ ) 3 ] | sin p ' | 

IK IK IK R 


IR 
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For U < 0: 


dF__ = [0.356 + 0. 321R - 0.368(R__) 2 + 0.392(R ) 3 ] 

RF IR IR IR 

(1 - |sinp^|) + [0.356 + 0.013lR ffi - 0.0764{R IR ) 2 (A-46) 

- 0.008 5(R ffi ) 3 ] |sinp^| 

dF pR = [0.356 - 0. 151R IF - 0.314{R if ) 2 + 0. 164(R if ) 3 ] (1 - |sin0j.|) 

+ [0.356 - 0. 0131R - 0. 0764(R ) 2 (A-47) 

1±* lb 

- 0.0085(R IF ) |sin0j,| 


where for each rotor 

Rj = tan 1 

r 

i x_ x i 

P' = tan * 

v 

H 




(A-48) 

(A-49) 


The thrust coefficients in Equation A-42, C TF and C TR , can be calculated 
through each rotor trim state, that is. 



(A-50) 


Although T can be calculated by trimming the previously derived equations, 
c 

Reference 2 has performed this analytically: 


T 

c 


(X - Xj) 
2 



+ IJ. \n 


0 

oo 


2 



0SO . 
~2~ i 


(A-51) 


The above equations represent a trancendental set for the solution of X __ 

lb 

and \ and involve both front and rear rotor trim states. 

IR 


99 



Trim Calculations 


Trim can be calculated by solving the previously derived equations, assuming 
the rigid body trim states are known. Reference 2 develops these equations 
for each rotor. 


Coning angle: 


8 


oo 12 


e u p 

4 T + °° + 1 - ■ °° 
c 6 5 2 


(A- 52) 


Longitudinal flap angle: 
4 


8 


co 


2-i 

1 - £— 

2 




(X 'V + 2 * T 3 . 

2 + T e oo + — ■ T ""so* 


SO 


Lateral flap angle: 

- 4 tL 


8 


so 


’ 1 + V 


2 -i^ oo + ^co 


(A- 53) 


{A- 54) 


Control Input Coupling 


The actual control inputs are collective (6 _), differential collective (8 ), 

O JLJ 

longitudinal cyclic (6 ), lateral cyclic (6 ) and differential cyclic (6 ), The 

Lj b JK 

derived control inputs 8, 0 , and 0 are related to these as follows; 

c s 


0^ = TDC + TDBfi^, 

OF C B 

(A-55) 

0 = TDC 6 - TDB8„ 

OR C B 

(A- 56) 
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e cF = TDH 6 r +TDS6 g 


(A- 57) 


e CR = TDR *R - TDS6 S 

Coupling to Rigid Body Equations 


(A-58) 


Although a complete set of equations of motion can proceed to include the full 

rigid body aerodynamics, the model was modified to utilize previously derived 

2 

rigid body equations. This allowed the use of trim data for the rigid body 

states, u , v , w , 6 , and 0 and the control inputs 5_ , 8_ , 6_ , and 
o o o o o Bo Co So 

S Ro* 

Some modifications to the force and moment equations are also required. 
Because the HELCOP model in Reference 2 consists of perturbations of rigid 
body states to rigid body equations, the rigid body perturbations of Equations 
A-28 through A-37 need not be performed. 

Deresidualization for Total Vehicle/Rotor Dynamics 


Because the Reference 2 data includes rotor trim states and rigid body 
dynamics coupled through the rotor, the required vehicle /rotor coupling 
involves proper treatment of rigid body dynamics which ’’pass through" the 
rotors. This procedure is best explained using a state space formulation. 

NASA has derived its HELCOP model in the following form: 

x = Fx + Gu (A-59) 

a a 
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where 


x is the rigid body state vector (eighth order) 

cl 


x = (u,v,w,p,q,r,9,0) 
a vv v v vvvv 


u is the control vector 


rp 


u 




F is the state coupler matrix 
G is the control coupler matrix 


The current development is represented as 



where: 


(A-60) 


x is the rotor state vector (12th order) 
r 


x 


0 


'OF* 0 OF' 8 


,B .B , B , B ,8 , B 

CF P CF P SF P SF P OR OR P CR' 


CR* 8 SR’ ^SR* 


Because HELCOP uses rotor trim states, the NASA model A- 59 is an implicit 
residualization of Equation A-60. Residualization of A-60 is performed by 
setting x^ to zero and solving for x^: 


x 

r 


= -F" 1 (F x +G u) 
r,r r, a a r 


(A-61) 
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This results in 


x 

a 


(F 


a, a 


F F _1 F )x +(G 
a, r r,r r, a a a 


F F" 1 G ) u 
a»r r,r r 


Comparing to Equation A- 15 yeilds 

F = F - F F" 1 F 

a, a a,r r,r r, a 


G = G - F F G 
a a, r r* r r 


(A-62) 

(A-63) 


The current situation, however, is peculiar in that we know F, F , F , 

a # r r, r 

F , G and G but not F . The solution for F can then be obtained 
r, a a, r a, a a, a 

using Equation A- 60 by "deresidulizing** F: 


F =F+F F -1 F 
a, a a, r r,r r,a 


(A- 64) 


This then provides the missing part of our formulation. Equation A- 16 
could also be used to specify G & ; but this is already calculated. The 
equation therefore provides a useful validity check on the rotor modeling 
task. 


Comparing G of Equation A-63 with the calculated control effectiveness 
values of HELCOP was performed for key input parameters. Figures A-4 
to A- 12 contain some comparisons. Minor discrepancies can be attributed 
to numerical accuracy. Errors at high forward velocity (particularly x = 

160 kn) are caused by the inability of the rotor modeling computer program 
to achieve perturbations on rotor interference parameters (see Equations 
B-42 and B-43 of Reference 2). These calculations are transcendental, and 
convergence on the wrong interference parameter for the 160- kn case results 
in gross errors in effectiveness values. For this reason, the rotor model 
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Figure A- 12. Comparisons 


is considered suspect at airspeeds higher than 120 kn. This should pose 
little problem in control design as sufficient data exists for gain scheduling 
in the remaining flight conditions. 

One interesting result concerns the value of L 6 in Figure A- 9. A key roll 

parameter, the HELCOP derived values were consistently lower than reference 

data (although not significantly). The residualized rotor model results were 

consistently higher, as shown. The sensitivity of L gg to the flapping spring 

constant, K ft , shows a definite impact of Lj . Since HELCOP does not 
p s 

include this in its calculation, one can conjecture that the cause of the 
discrepancy between calculated and referenced results could be the effect 
of the flapping spring constant. One could further conjecture that the value 
of K = 10000 (an estimate used here) is too high. 

p 
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One important parameter not shown is 75 • Because of the scaling used in 

s 

Figure 13 of Reference 2, the current results could not be plotted. However, 
all values were approximately 15% low. 

Rotor Measurements for Feedback Control 


Selected designs in Section 3 utilize states of the rotor models for coning, 
P_._ and p__, and lateral flapping, P e _ and P__, to increase the control 
bandwidth. Because these are harmonic coefficients, they must be derived 
from the appropriate geometry of actual measurements. If we have a sensor 
measuring total flap angle on each blade of a three-bladed rotor, the appro- 
priate states can be derived. 

Let; 

Pi = " e c cos ♦ ” e s Sin ’* f 

p = P n - P cos (t + 120°) - 0 sin (i|f + 120°) 

&L u c s 

p„ = p n - P cos (f - 120°) - B sin (* - 120°) 

o U C S 


The simple algebra and trigonometry results in 


\ 

TD 

1 


cost. -sin 


" B 2 + P 3' 2S l" 


1] 



CO 



sint, cost 


T (P 3 " V 


This also requires that t be measured. 


This, of course, ignores higher harmonics, blade lagging, and flexure. 
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APPENDIX B 


MODAL CONTROL DESIGN TECHNIQUES 


Introduction 


Modal control design is one of placing closed loop eigenvalues and eigen- 
vectors for a given plant. The term modal refers to the mode design (pole 
placement) and the construction of system responses to these modes. The 
result is a design in which various system responses contain dynamics 
associated with desired roots or modes. For example, the lateral- directional 
axis of a fixed wing aircraft has three dominant modes, roll, dutch roll, and 
spiral. Pole placement algorithms can be used to change the roots of the 
closed loop system corresponding to each of these modes; however, unless 
aircraft responses are isolated from certain modes, the transient response 
of the system might not demonstrate the desired properties. The dutch roll 
root should dominate the (3 - R responses, the roll root should dominate the 
P response, and the spiral root should dominate the 0 response. 

Most aircraft contain some natural response decoupling, but existing pole 
placement algorithms do not account for this. System eigenvector placement 
provides the key additional design element. 

Modal Control Overview 


6 

The design theory, based upon work done by Moore, may be presented 
briefly as follows. Using state space notation 
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x = Fx + Gu 


y = Cx 


(B-l) 


where 

x(t)eR n ; u(t)eR m ; y(t)eR P 
F, G, C have the appropriate dimensions and 
rank (G) = m 
rank (C) = p ^ n 


An eigenvalue- eigenvector pair for a closed loop design for system B-l is 
defined as 

F v. = X.v. (B-2 ) 

c x IX 

where X^ is one of the n eigenvalues, v. is the corresponding eigenvector, 
and F c is the closed loop plant matrix. 


The X., v. pair can be achieved with control through the input matrix G; 
that is, 

Fv. + Gw. = X.v. (B-3 ) 

x x x x 

where w^ is a vector which satisfies the equation for the desired X^ and v^. 


Case I: State Feedback 


If we have access to all the states of the system, the feedback control law 
becomes 
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u = Kx 

(assuming p=n and letting C=I without loss in generality) 
and B-2 becomes 

(F + GK)v. = X.v. (B-4) 

Matching B-4 with B-3, we must find w^ and K, which maps into w.. 

w. = Kv. 

1 1 

For the entire collection of closed loop eigenvalues and eigenvectors, A-3 
becomes 

VA - FV = GW (B-5) 


where 

A - diag (X^ X 2 , . . . X n ) (nxn) 

V is the eigenvector matrix (nxn) 

W is the matrix which satisifes B-5 for A and V (mxm). 


Likewise from B-4, 


W = KV 


(B-6) 


therefore 

K = WV' 1 * (mxn) 


(B-7) 


* 6 

Moore demonstrates that rotating complex eigenvector pairs to real vectors 
produces a real K matrix which performs the desired placement. 
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Case 2: Output Feedback (p<n) 


y = Cx 

In this case Equation B-6 becomes 
W = KCV 


(B-8) 


Since C is of rank p<n, a unique solution for K is impossible. The alternative 

* 

used here is to choose only p eigenvalue-eigenvectors to place; in other 
words, choose pX’s and V (nxp) and solve for K such that 

K = W(CVf 1 (mxp) (B-9) 


Design Constraints 


The key design issue is to find W, which satisifes B-8. In general, one 
cannot completely satisfy both exact eigenvalue and eigenvector placement. 

Case 1: Single Input Systems 

For single inputs Equation A-3 reduces to 
<¥ ' F) v i = 

where is a scalar. Since this single variable must be adjusted to place 
n parameters on the left hand side, X^ and n-1 parameters of v^, the job is 
impossible. Therefore, the single input case only involves pole placement 
with arbitrary eigenvector position. 


There is no guarantee that the remaining closed loop system eigenvalues 
will be stable. 
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Case 2: Multiple Input Systems 


The real benefit of modal control is employed when more than one control 
is available. The job becomes a matter of choosing the right w^'s in 
Equation A- 3. 

(X.I - F)v. = Gw. 

1 l l 

The design procedure initially involves choosing portions of v^ in order to 
eliminate certain state responses from a mode while emphasizing others 
and letting other responses (control or compensation) react arbitrarily. 


For Rank (B) = m, m free parameters can be specified, one of which is the 
eigenvalue. Equation B-3 can be rewritten 


F ci q t = t (F “ q i = 0 n b y < n+m ) 


V 

i 



(B-10) 


q^ is therefore a null space mapping of F^. A convenient tool for finding 
the relationship between v. and w. is contained in the singular value 

g 1 ^ 

decomposition of F^. 


The singular values of the matrix A . are the eigenvalues of [A .A .*], 

OX ox ox 

where * refers to the complex conjugate transpose. They are observed 
through the singular value decomposition of F c .. 


F . 
ci 


X.E.Z. 

ill 


(B-ll) 



where: 


X is an n-by-n matrix containing columns of orthogonal left singular 
1 

vectors of F . 

ci 

E,. is ann-by-n+m matrix containing n singular values, o's, of A . 



= [E^, [0]]; E^ is nxn diagonal 

Z., is an n+m-by-n+m matrix containing n+m orthogonal right singular 

vectors of F . 

ci 

By rearranging B-ll 

F .Z. = X.E. (B-12) 

ci 1 1 1 

and noting that the last m columns of the X.E^ product are null, we find the 

appropriate null space for F . by using the last m columns of Z.: 

OX X 

F . Z. = 0 (B-13) 

Cl 1 


115 



where Z. is defined as 
1 

Z. = [Z. Z.] ) n+m- (B-14) 

1 ^ ^ 


The matrix Z. is a set of m orthonormal basis vectors spanning the null 
space of A Referring to B-10, we have 


q i 


= Z.a. 
i i 


(B-15) 


where is an m vector of linear coefficients not all of which can be zero. 

l 

Modal Control Algorithm 


Returning to our original problem, we wish to specify the desired v. and 
find a which solves B-10. This cannot be done in general if m<n, there- 
fore we formulate a least squares performance index 

J = (v di " V Q(v di " v i* (B-16) 

where 

v^. is the desired eigenvector 

v. is the resulting admissible vector subject to the constraint B-10. 

Replacing v^ in B-15 with 
v. = E Z.a . 

i li 
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where 


E is an n-by-n+m matrix 



n m 


and minimizing B-16 with respect to at 


= Argfmin J ] 

.*•_» ^ r T «- 1 * nr 

= (Z. E Q EZ.) Z. E Q V di 


the appropriate W. is found using B-15. 


>;< x 


-1 _* T 


W. = E Z.a. = EZ.(Z, E Q EZ.) Z.E Q 

1 XI X X ^ 1 l 


V 


di 


where E is an m-by-n+m matrix 



n m 


(B-17) 


(B-18) 


Referring to B-16, it should be noted that one could try to fit the n elements 
of to a prespecified n elements of Unless happens to lie in the 

m-dimensional subspace of a perfect fit is impossible. As a practical 
matter one should specify up to m elements of and let the other elements 
be arbitrarily placed. This is performed by the proper choice of Q, that is, 
by placing the desired quadratic weights only in diagonal elements corre- 
sponding to element locations in which are to be placed. The example 
which follows demonstrates this. 
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Illustrative Example: F-4 Lateral- Directional Axis 


The design of an inner loop control law for the F-4 lateral- directional axis 

was undertaken to illustrate the multivariable optimal control design 

7 

procedures developed under an ONR sponsored contract. Much of the 
problem description is repeated here because the eigenvalue- eigenvector 
design goals are identical. 


The F-4 fighter aircraft at a low dynamic pressure flight condition is taken 
from Reference 8. The dynamics are 

x = Fx + Gu 


with 


x = 



stability axis roll rate 
stability axis yaw rate 

f 

angle of sideslip 
bank angle 
rudder deflection 
aileron deflection 


u = 


rc 


ac 


rudder, command 


aileron command 


118 



Matrices F and G are 


F = 


G 


-.746 

.387 

-12.9 0 . | .952 

6. 05 

.024 

-.174 

4.31 0. 1 -1.76 

I 

-.416 

.006 

-.9994 

-.0578 .0369 .0092 

-.0012 

1. 

0. 

• 

o 

• 

o 

* 

o 

0. 

0. 

0. 

1 

0. 01 . j - 20 . 

0. 

0. 

0. 

O 

• 

o 

* 

o 

• 

- 10 . 

0. 

0. 



0. 

0. 

Open Loop Poles 


0. 

0. 

X roll subsidence 

= -.079 

0 . 

0. 

X dutch roll 

= -.098 +j2. 

20 . 

0. 

X spiral 

= -.0063 

0. 

10 . 

X rudder actuator 

= - 20.0 


_ 

X aileron actuator 

* — 1 0 * 0 


An initial condition response for the open loop system is shown in Figure B-l. 


From the point of view of fighter handling qualities, all four of the lateral 
axis closed loop roots have desired values which can be taken from MIL- 
F8785B, as is done, for example, in Reference 9. The desired roots are: 

a) Roll subsidence mode = -4. 0 

b) Dutch roll mode = -0.63+j2.42 

c) Spiral mode =-0.05 
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P (rad/sec) 
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0 1 2 3 4 5 

Time (sec) 

Figure B-l. F-4 Lateral- Directional Open Loop Response 


Each of these poles can be assigned an asymptotic eigenvector, v., which 

distributes the modal response, eV, among the state variables and outputs 

of the system. However, each eigenvector is constrained to lie in a two- 

dimensional subspace (page 16, Reference 7). An element of this subspace 

was selected by using B-18 to find the best linear projection of an unconstrained 

desired vector, v,., on the subspace. The results of the eigenvector 
di 

selection are: 

-4t 

a) Roll subsidence mode (e v^ 

Desired v„ = [1. 0 0 a a a] 

dl L 

Attainable v^ = [1. -.007 0 -.25 . 13 -. 56] 
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0 (rad) 




**• S3 1 

b) Dutch roll mode, real part e ’ (cos2.42t)v 2 

Desired v^g = [0 a 1, 0 a a] 

Attainable v 2 = [0 15.6 1. 0 7.86 -.103] 

f -.631 

c) Dutch roll mode, imaginary part e * (sin 2.42t)Vg 

Desired = [0 1, a 0 a a] 

Attainable Vg = [0 1. 6.16 0 -9.49 14.6] 

d) Spiral mode | e *^^v4 

Desired v^ = [ a a 0 1. a a] 

Attainable = [-.05 .037 0 1. -.0014 -.0079] 

A few comments are in order to explain these choices. Consider, for 
example, the roll subsidence mode. The desired eigenvector is taken to be 
= (1 0 0 a a a), which means that the mode should show up dominantly on 
roll rate, but not on yaw rate or sideslip (we want no sideslip buildup during 
turn entries). These are good basic handling quality considerations. The 
a's in the vector indicate that we do not care how much of the mode shows up 
on these components. Certainly, since 0 ~ Jp dt, some mode content has to 
be expected on element a^ and, similarly, if the surfaces are actually con- 
trolling the mode, some mode content should also appear in a_ and a . The 

5 6 

linear projection which best achieves these objectives is shown as v^ above. 
Note that we can satisfy our desires almost perfectly. 

Similar arguments also apply to the dutch roll mode. In this case we want 

no oscillatory dutch roll content on roll rate and bank angle. This is a key 

8 

handling quality requirement for all well-behaved lateral control laws. 
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This output feedback problem poses little difficulty using modal control 
techniques since we wish only to place four eigenvalue- eigenvector sets. 
Actually there is concern for the remaining two poles. It is desirable to 
avoid moving these remaining poles too far to the left, which would hint at 
high actuator demands. The control design gives us a bonus here if the 
measurement vector contains no actuator states; in other words. 


y 


[C, 0] 


X vehicle 
X actuator 


Since the form of the open loop system matrix F is 

F , F 
11 12 


0 , F 


22 


then the closed loop system matrix F £ is 


F . F 
11 ' 12 


F = 
c 


BKC, F 


22 


The trace of a matrix is the sum of the eigenvalues; therefore. 


n n 

Y X . = trace F = T* X . = trace F 
L, ci c i 

i=l 1=1 


The sum of the closed loop eigenvalues is equal to the sum of the open loop 
eigenvalues. This becomes important if our design goal involves moving 
plant poles to the left (which is usually the case), because it requires that 
the remaining eigenvalues of the closed loop system move to the right. Since 
the "remaining" eigenvalues are our "actuator poles" (roughly speaking), we 
are assured that they do not increase in magnitude. 
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Design Results 


Case I: Vehicle State Feedback 

Using the vehicle states listed earlier, 

T 

y = (p s » r s » 0. 0) 

Perfect pole placement was achieved (as guaranteed by earlier analysis) 
and the eigenvectors which we found to be achievable were obtained. The 
remaining poles move to the right as expected. 

X 5 = -19.03 (formerly -20) 

X„ = - 6.64 (formerly -10) 

D 

An initial condition response, shown in Figure B-2, verifies the eigenvector 

response isolation desired. This is virtually identical to the results of the 
7 

ONR design study, which used optimal feedback control, although some 
minor pole movement was observed in the later design when actuator gains 
were eliminated. 


Case II: Practical Output Feedback 

Although 0 feedback is possible it is usually not advisable because of the 
difficulty of obtaining good side slip data. A more reasonable approach is 
to use lateral acceleration instead. This is not a trivial linear transformation, 
since the measurement includes actuator states. Also, our nice property 
of invariant eigenvalue sums is no longer valid. 
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Using the output vector 


T 

y = (p , r , n , 0) 

J s s y 

we achieved the same desired "vehicle" poles and eigenvectors as in Case I 
(as we should) plus two additional poles at 

X c = -15.82 

D 

\ 6 = - 5.98 

Although not guaranteed, the two additional poles move even further to the 
right than in Case I. 

The initial condition response is virtually identical to that in Case I (Figure 
B-2). 


Case III: Reduced Measurement Set 

Another practical alternative is to eliminate more sensors (in other words, 
save $'s) and go for the same design goals. One logical candidate would be 
the roll attitude gyro (measuring 0). In modal control this would require 
placing one fewer pole (and eigenvector). The logical candidate here would 
be the spiral mode (X spiral desired = - .05). 

Results were impressive, with the initial condition response (Figure B-3) 
showing little impact from this change. A slight decrease in dutch roll 
damping is observed when Figures B-2 and B-3 are overlaid. Additional 
eigenvalues are 
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2 


Time (sec) 


3 


4 


5 


B-3. Closed Loop Design: Case 3 




023 


^ spiral 

X, =-15.83 

5 

X c = - 5.99 

b 

The spiral root does not meet our objective of - .05; however, it is unlikely 
that this value would produce significant comments from rating pilots. Even 
slightly unstable spiral roots are usually acceptable. 

Further Measurement Reductions 


If our sensor coffer is poor and we need to eliminate yet another measure- 
ment (or the last one of a kind fails) let us see what we get. 

Case IV: y T = (r, n ) 

y 

In this case the logical design goal to eliminate is the roll subsidence mode, 

since we are no longer using the roll rate gyro to measure p . 

s 


Design results place the dutch roll pair exactly and produce the following 
remaining poles: 


^roll subsidence 

^spiral 

X 5 

X 6 


- 2.54 (-4.0 was our earlier goal) 

- .296 (-.05 was our earlier goal) 
-16.00 

- 7.27 
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The transient response is shown in Figure B-4, The impact of the roll 
subsidence and spiral modes on the dutch roll is obvious (decreased damping). 
The important factors are: 

1, The dutch roll oscillation is not contained in the roll 
response because the dutch roll eigenvector placements 
preclude this. 

2. The response is still a vast improvement over the open loop case. 

Case V: y T = (p g , r g ) 

This case is much more difficult because no logical fallback on our design 
goals exist, that is, we cannot place half of a complex pair. Therefore, it 
was decided to try and place the dutch roll mode (as in Case IV). 

Modal control guarantees the prescribed pole and eigenvector placement 
(since our earlier analysis showed it was feasible). Other results, however, 
demonstrate the disadvantages of this: 

X 3 4 = - 4. 62 +j 19, 63 

X, = - .0047 

5 

X c = -20.47 

b 

The transient response is not shown because the design is obviously 
unacceptable. The point here is that an unreasonable task was required, 
given the current sensor /hardware configuration. (We actually got what we 
asked for, but should we have asked?) 
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Conclusions 


In this appendix modal control is developed from its conceptual foundation 
(provided in Reference 6) to a practical design tool through partial eigen- 
vector placement techniques (developed in Reference 7). The F-4 lateral- 
directional axis design example demonstrates the usefulness of this design 
method to perform precision mode shaping and tradeoff studies. 
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A PPENDIX C 


CH-47 MODEL TESTS 


Introduction 


This appendix contains results of the control systems tests conducted on the 
NASA-Langley research vehicle CH-47B. Included are static and dynamic 
tests conducted in the hangar plus results of tests conducted in hover and 
cruise flight. 

The purpose of these tests was to provide a check on the dynamics of various 
control system elements, including the presence of any nonlinearities. A 
reasonably accurate model of the control system is needed to permit continu- 
ing investigation into ways and means of improving the gain and bandwidth 
capabilities of the system. 

Static and dynamic measurements on some portions of the system could be 
conducted in the hangar with hydraulic and electric power-on but rotors 
stationary. A function generator (sine, square, or triangle) was used to 
introduce signals of appropriate amplitude and frequency into the ECS servo. 

A seven-channel strip chart recorder was used to record ECS servo position 
and/or upper boost actuator position, as well as the forcing function. For 
the hysteresis tests, an X-Y recorder was used. 

It was not possible to measure the frequency response to include the upper 
boost actuators in the hangar because of the excessive flopping of the rotor 
blades produced in the non- rotating, unloaded condition. These were obtained 
under hover and cruise flight conditions, with the results recorded by the 
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special on-board instrumentation and magnetic tape recording system. 

Static control system perturbation gains were obtained from the low frequency 
response data at hover and corrected for hysteresis. 

During all tests, the normal ships SAS system was disabled to prevent inter- 
ference with the measurements. This was accomplished by making changes 
to the SAS servo wiring so that one servo was caused to go hardover in the 
fully extended position while the other servo was fully retracted. They 
therefore cancelled one another in the linkage while each still acted as a 
rigid link, with no affect on the control system dynamics. 

CH-47B Control Characteristics 

Nominal Control Characteristics 

The CH-47 flight control system uses mechanical linkages, hydraulic boost 
actuators, and complex mechanical mixers to transmit pilot control motions 
to the rotor swash plates. The resulting control system model can be repre- 
sented by a relatively simple block diagram, as shown in Figure C-l. This 
diagram includes effects of backlash in the linkage and presents the static 
and dynamic relationships and scale factors determined in this investigation. 

The ECS servo is not part of a normal CH-47 control system but was added 
to this research vehicle in a prior program. It permits the introduction of 
electrical commands from whatever source desired into the control system. 

It operates as a parallel servo with the safety pilot's manual controls when 
a connecting clutch is engaged. This enables the safety pilot to quickly 
regain control by de- clutching the ECS servo when conditions warrant. 
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The gains shown relate ECS or safety pilot control stick motion to the net 
sum of the upper boost actuator motions by axis. For example, 1 cm roll 
ECS servo motion produces 3. 9 cm of delta roll; this is the summation, as 
is shown in Figure C-l of the four upper boost actuators in response to roll 
commands. 


ECS Servo Characteristics 

The ECS servos in the pitch, roll, and yaw axes were each checked at three 
amplitudes for frequency response, rate limits, and nonlinearities. The ECS 
was de- clutched from the rest of the control system during these tests. 

Measured Frequency Responses. -Frequency response data was recorded 
for the range 0. 25 to 20 Hz and at peak amplitudes of +2. 5%, +5%, and +10% 
of the servo stroke of +5. 08 cm (+2 inches). Rate limiting was encountered 
at the higher frequencies, depending upon stroke amplitude. The results 
of these measurements are presented in Table C-l. 

No Linearities. - Linearity and backlash of the ECS servo was evaluated 
by driving it with a 0. 1 Hertz triangle wave and recording input vs. output 
on an X-Y recorder. The results are presented in Figures C-2 and C-3. It 
should be noted that the transfer is essentially a straight line with negligible 
gap between the forward and return stroke. 

Rate Limits. -During the frequency response tests, rate limiting was 
encountered at the 5% and 10% strokes at the higher frequencies. Analysis 
of these conditions revealed the rate limits shown in Table C-2. 
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TABLE C-l. ECS SERVO FREQUENCY RESPONSE DATA 


INPUT 

PITCH 

ROI.I. 

YAW 

Frequency 

Amplitude 

Amplitude 

Phase 

Amplitude 

Phase 

Amplitude 

Phase 

Hertz 

% Stroke 

DB 

Deg. 

DB 

Deg. 

DB 

Deg. 

. 25 

2. 5 

0 

0 

0 

0 

0 

0 

.5 


0 

0 

0 

0 

0 

0 

1 


0 

0 

0 

0 

0 

0 

2. 5 


2 

-10 

-. 2 

-5 

-.2 

-10 

D • 


2 

-20 

-. 28 . 

-15 

-. 36 

-20 

10. 


-1. 1 

-35 

-. 98 

O 

CO 

t 

-. 64 

-25 

2 0. 


-4.4 

-55 

-3. 6 

-55 

-3.4 

-50 

. 25 

5 

0 

0 

0 

0 

0 

0 

. 5 


0 

0 

0 

0 

0 

0 

1. 


0 

0 

0 

0 

0 

0 

2. 5 


0 

-5 

-.2 

-5 

-. 1 

-5 

5. 


0 

” 1 0 

-.2 

- 1 5 

-. 18 

-1 5 

10. 


-3. 3 

-55 

-3. 0 

-55 

-2. 9 

- 55 



rat 

e limit 

rat 

e limit 

rat 

e limit 

2.5 

10 

0 

0 

0 

0 

0 

0 

. 5 


0 

0 

0 

0 

0 

0 

1. 


0 

0 

0 

0 

0 

0 

2. 5 


0 

-5 

- 9 

♦ 

- 5 

0 

-5 

. 5 


-3. 1 

-45 

-2. 7 

-45 

-3 

-45 



rat 

1 

e limit 

1 

rat 

0 limit 

1 

rat 

1 

e limit 

i 


Notes: Amplitude ratios and phase read from strip chart recordings. Phase stated 
to nearest 5 deg. 
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TABLE 02. ECS RATE LIMITS 


Amplitude 

Rate Limit, cm/ sec (in/ sec) 

% 

Pitch 

Roll 

Yaw 

2. 5 

7. 62 (3. 0) 

7. 62 (3. 0) 

7. 62 (3. 0) 

5 

7.37 (2.9) 

7. 62 (3. 0) 

7. 62 (3. 0) 

10 

7.11 (2.8) 

7.87 (3.1) 

7.87 (3.1) 


Math Model. -From the foregoing results we conclude that the ECS 
servo can be represented by a first-order lag with a corner frequency of 
approximately 15 Hz. This appears to be valid up to a frequency of 20 Hz, 
the highest frequency investigated. As an illustration, the data from Table 
C-2 for the pitch servo at 2, 5% stroke is plotted in Figure C-3. Also 
included are the phase and amplitude curves for a true first-order lag with 
a corner frequency of 15 Hz. 

Primary Control System Characteristics 

The primary control system evaluated here includes that portion from the 
output of the ECS servos to the output of the upper boost actuators, which 
drive the swash plate. Included between these limits are the lower boost 
actuator and the mechanical linkages and mixers. 

It was not possible to conduct frequency response tests on the primary 
control system in the hangar. This was because the oscillatory motion of 
the upper boost actuators tended to start the stationary rotor blades bouncing, 
creating a condition possibly hazardous to blades or personnel. Therefore, 
frequency response data was only obtainable during flight. However, it 
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was possible to obtain some data, including nonlinearities, step responses, 
and rate limits, in the hangar. 

Measured Frequency Response . -Frequency response data was recorded 
in flight on magnetic tape by the on-board instrumentation system. The 
frequency range investigated was 0.25 to 8 Hz and the amplitude +10% of ECS 
stroke, except at the higher frequencies of 4 and 8 Hz, where amplitudes of 
+ 5% and +2. 5% respectively were used to avoid rate limiting in the ECS 
actuator. The four upper boost actuators were added algebraically according 
to the formulas given in Figure C-l. This gave the net pitching, rolling, or 
yawing output from the control system. The digitized data was then subjected 
to a Fourier analysis program which was able to provide the amplitude and 
phase shift of the funamental forcing frequency as well as the noise frequency 
components in the sensor outputs. 

Fourier analysis was then used to determine the transfer function of the 
control system from ECS servo output to the summation of the upper boost 
actuators on a per-axis basis. The results of this analysis are presented 
in Table C-3. 

Nonlinearities . - Linearity and hysteresis of the primary control system 
were checked by driving the ECS servo with a 0.1 Hertz triangle wave and 
recording the ECS servo output and the summation of the upper boost actuator 
outputs on an X-Y recorder. The proper summations are listed in Figure 
C-l. These tests were conducted at three ECS servo output amplitudes, 
+2.5%, +5%, and +10% of full stroke. The results for the pitch, roll, and 
yaw axes, respectively, are shown in Figures C-4, C-5, and C-6. 
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±5% STROKE 



Figure C-4. Pitch Axis Hysteresis 
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Figure C-6. Yaw Axis Hysteresis 


. 78 % 


The figures show that there is considerable hysteresis in all three axes, 
with pitch exhibiting the most. The magnitudes of the hysteresis as a per- 
centage of full ECS servo stroke have been extracted from the figures and 
listed in Table C-4. 

It should be noted that, at least over the range of control system motion 
induced by ECS servo strokes up to +10%, the system is quite linear, that 
is, linear in the sense that the graphs have very little curvature. They 
were also quite repeatable since several traces of the hysteresis loop fell 
on top of each other. 

To test for control system linearity at positions other than control stick 
dead center, the roll stick was moved off center about 50% of its travel and 
the roll hysteresis loop repeated. The results were very close to the results 
obtained for a zero centered stick. 

To check for frequency- dependent effects in the hysteresis loops (for 
example, spring- dashpot effect), the roll axis tests were repeated using 
other velocities for the forcing function triangle wave. The frequency of the 

TABLE C-4. HYSTERESIS (% of full travel) 


ECS Stroke 

Pitch 

Roll 

Yaw 

+ 2.5% 

1.1 

.54 

. 72 

+ 5.0% 

1.3 

.64 

.78 

+ 10% 

1.4 

. 62 

.89 
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wave was varied from 1/5 Hz to 1/30 Hz tests. It is concluded, therefore, 
that the hysteresis is due primarily to simple mechanical backlash or 
compliance plus friction. 

Step Response. - Since no frequency response tests of the control system 
were possible in the hangar, as explained earlier, a step response measure- 
ment of the pitch axis was substituted. A manual transient input was applied 
to the pitch stick, using an improvised solid stop to control the amplitude. 

The summation yielding pitch was recorded on the strip chart recorder. The 
result of this measurement is replotted as the solid curve in Figure C-7. 

The dashed curve is a plot of a pure second- order lag with a natural frequency 
of 50 rps and a damping ratio of 0. 65. It should be noted that this is very 
close to the transfer function of the upper boost actuators given in Reference 1. 
This was 50 rps and damping of 0. 55. It also approximates the third-order 
model obtained by frequency response testing discussed later (and included 
in the pitch block diagram of Figure C-l). 

Rate Limits. -The rate limits of the combined lower boost and upper 
boost actuators were determined by applying step inputs to the control stick 
and yaw pedals, one axis at a time, and measuring the resulting upper boost 
actuator position with time on the strip chart recorder. From the slope of 
the curve, the maximum rate was calculated. While that was not considered 
to be a sophisticated measurement, the purpose was to determine the lower 
and upper boost actuator combination rate compared with the previously 
determined value of 7. 62 cm/s (3 in/s) for the ECS servo. This would 
determine which element was controlling in terms of rate limiting. The 
results of this test are tabulated in Table C-5. 
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TABLE C-5. UPPER BOOST ACTUATOR RATE LIMITS 



Rate Limit, cm/s (in/s) 


Pitch 

Roll 

Yaw 

At boost output 

23.4 (9.2) 

71.4 (28.1) 

57. 9 (22.8) 

Equiv. ECS axis 

14.0 (5.5) 

18.8 ( 7.4) 

11.7 ( 4.6) 


In each case the equivalent ECS servo rate limit of the lower and upper boost 
actuator combination exceeded the actual ECS rate limit of 7.62 cm/s (3 in/s). 
Therefore, the ECS servo is the rate limiting element in the system. 

Math Model. -From the results of the hysteresis and frequency response 
tests, a model for the control system was derived in terms of idealized 
elements. 

The hysteresis values determined from the 10% amplitude tests were applied 
to the 5% and 2.5% amplitudes as well, since the best measurement was 
obtained at the largest amplitude. Also, from other tests it was known that 
a second-order lag at approximately 50 rad/sec with damping about 0. 6 
could be expected. The results of synthesizing the control system with these 
elements are presented in Table C-6. It was necessary to include a first- 
order lag at 100 rad/sec to provide a reasonably good match of the recorded 
data. This must represent the dynamics of the lower boost actuator. In 
order to provide the best curve fit, the natural frequency and damping ratio 
of the second- order lag were varied somewhat from axis to axis. No 
explanation is available for this discrepancy, since the same upper boost 
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TABLE C-6. CH-47B CONTROL SYSTEM RESPONSE CURVE FITTING 
(Hover Data) 


Hysteresis 


Second Order 


First Order 


Gain 

Phase 

Deg. 

» r 

1 , 4% of 

4 inches 

.93 

-10 

.93 

-10 

.93 

-10 

.83 

-20 

.47 

-40 

.62% oi 

4 inches 

.98 

- 4 

.98 

- 4 

.98 

- 4 

.98 

- 4 

.98 

- 4 

.85 

-16 

. 89% Of 

4 inches 


Phase 

Deg. 


Gain 


.45 Wc = 


Phase 

Deg. 


100 rps 

- 2 

- 4 

- 8 
-15 
-27 


Gain | 

Relative 

Normalized 

.93 

1 

.94 

1.01 

.96 

1.03 

.94 

1.01 

.41 

.44 

.98 

1 

.98 

1 

.99 

1,01 

1.0 

1.02 

1.06 

1.08 

.75 

.77 

.96 

1 

.96 

1 

,97 

1.01 

1.0 

1.04 

1.01 

1.06 

.43 

.44 


Flight Test 
Measurement 


Phase Normalized phase 
Deg. Gain Deg. 



























actuators enter into the composite axis displacement. If it had been possible 
to use larger amplitude displacement at the ECS without encountering rate 
limiting at the higher frequencies, the quality of the data might have been 
improved. To average the results obtained, it could be said that a second- 
order lag at 45 rad/sec and 0. 5 damping, plus a first-order lag at 100 rad/ 
sec represent the control system between ECS output and total upper boost 
actuator displacement, excluding hysteresis. 

Airframe/Rotor Characteristics 


From the recordings of pitch, roll, and yaw rates plus lateral acceleration, 
a measure of the airframe transfer function can be obtained. With the 
Fourier analysis also applied to these sensor outputs, an added bonus is the 
ability to determine the relative amplitudes of noise components at the con- 
trol frequencies of interest. 

Measure Frequency Responses 

From the Fourier analysis of the sensor output data, transfer functions were 
obtained relating the magnitude and phase with respect to the axis summation 
of the upper boost actuators. This data is presented in Table C-7. Plots of 
the data are presented in Figures C-8 through C-ll. In addition, these plots 
include the amplitude and phase of the upper boost actuator summations with 
respect to their ECS servos, taken from Table C-3. It should be noted that 
there is good agreement between hover and cruise conditions for the upper 
boost summations with respect to their ECS and also for the body rate sensors 
up to 2 Hz. Except for the pitch axis, the body rate data above 2 Hz appear 
to be unreliable. This is undoubtedly due to the high noise content of the 
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TABLE 07. AIRFRAME FREQUENCY RESPONSE DATA 


Flight 

Condition 

Forcing 

Freq. 

Hertz 

Delta 

4 

Pitch 

p 

Delta Roll 

r 

Delta Roll 

Ny 

Delta Roll 

P 

Delta Yaw 

r 

Delta Yaw 

Ny 

Delta Yaw 

Gain 

Phase 

Gain 

Phase 

Coin 

Ph as e 

Gain 

Phase 

Gain 

Phase 

Gain 

Phase 

Gain 

Phase 

Hover 

.25 

NR 


.138 

- 46 

.019 

- 35 

.138 

- 45 

.009 

-232 

.033 

- 91 

.0089 

-14 



.5 

.18 

- 76 

.093 

- 82 

.013 

- 88 

. 027 

-154 

.006 

-255 

.017 

-124 

.011 

+ 2 



1 

.089 

- 94 

.042 

-106 

.005 

- 99 

.016 

-198 

. 003 

-269 

.006 

-110 

.0076 

-12 



2 

.083 

-112 

.0076 

-105 

. 0018 

-245 

.0114 

-198 

.0054 

-223 

.0023 

- 79 

.0085 

+34 

y 

r 

4 

.04 

-199 

.032 

-151 

.0115 

-450 

.049 

-308 

o 

H-* 

CO 

-300 

.015 

-124 

.15 

-21 

41. 1 

8 m/s 

8 

.007 

-317 

.017 

-267 

.0038 

-483 

.134 

-535 

.0093 

-547 

,0005 

-106 

.012 

- 9 

(80 kn) 

.25 

NR 


. 137 

- 72 

.02 

- 62 

.035 

-122 

.035 

-132 

.028 

-115 

.0061 


Cruise 

. 5 

.141 

- 68 

.076 

- 79 

.0075 

- 53 

.023 

-168 

.01 

-280 

.013 

-104 

.0084 

0 



1 

.095 

- 94 

.038 

-109 

.0033 

-107 

.02 

-108 

.0037 

-325 

.0055 

-107 

.0099 

-14 

'I 

f 

2 

.064 

- ! 14 

.012 

-157 

.0019 

-175 

.011 

-230 

.0014 

-613 

.0029 

-160 

.0135 

-39 


Gains in radians /sec per inch or g's per inch 
Phase in degrees 
NR --not run 
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sensor outputs at the rotor frequency of 3.75 Hz and multiples thereof, which 
are close to and influence the data at the forcing frequencies of 4 and 8 Hz. 
This was not the case in the pitch axis where the pitch rate gyro data follow 
a more natural trend at 4 and 8 Hz as shown in Figure C-8. Examination of 
the spectral magnitude plots of the pitch rate gyro output, to be described 
later, show that the rotor frequency noise is less intense than in the other 
axes sensors. 

It is apparent from the response plots that we would probably have obtained 
much better data at the higher frequencies if we had chosen forcing frequencies 
farther away from the known rotor frequencies; in other words, the 4 and 8 
Hz values were poor choices. 

Sensor Noise Components 

The Fourier analysis of the sensor outputs to determine gains and phase 
shift at the test frequencies also yields all other frequency components present 
up to and the cut-off frequency used in digitizing the data. This provided an 
excellent opportunity to assess the absolute magnitudes of the rotor frequency 
components picked up by the sensors. It has been recognized by various 
investigators that sensor noise has probably been the limiting factor in 
obtaining higher gain control loops on this helicopter. Spectrum plots are 
shown in Figures C-12 through C-15. These all represent the case for a 
0. 5 Hz forcing frequency into the roll axis, hover condition. 
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Roll Rate Spectrum Magnitude (rad/$ec) 






Lat Accel Spectrum Magnitude (G'$) 



The four figures represent the pitch, roll, and yaw rate gyro outputs and 
the lateral accelerometer. It should be noted that the pitch rate is the 
cleanest of the four signals and the roll rate is probably the noisiest. Also 
in some cases, the peaks at frequencies of 3 per rev and 6 per rev exceed 
those at the fundamental rotor frequency (1 per rev). 
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